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PREFACE 


The  work  described  in  this  report  was  conducted  for  the  US  Army  Engineer 
District,  St.  Louis,  by  the  Concrete  Technology  Division  (CTD)  of  the  Struc¬ 
tures  Laboratory  (SL)  and  the  Information  Research  Division  (IRD)  of  the 
Information  Technology  Laboratory  (ITL) ,  US  Army  Engineer  Waterways  Experiment 
Station  (WES).  The  investigation  was  authorized  by  DA  Form  2544,  Intra-Army 
Order  for  Reimbursible  Services,  No.  ED84-12,  dated  December  1983. 

The  investigation  was  accomplished  under  the  general  supervision  of 
Messrs.  Bryant  Mather,  Chief,  SL;  James  T.  Ballard,  Assistant  Chief,  SL; 

John  M.  Scanlon,  Chief,  CTD;  and  Dr.  N.  Radhakrishnan ,  Chief,  ITL,  and 
under  the  direct  supervision  of  Mr.  C.  Dean  Norman,  Program  Manager.  This 
report  was  written  by  Messrs.  Anthony  A.  Bombich  and  C.  Dean  Norman,  CTD,  and 
H.  Wayne  Jones,  ITL. 

COL  Dwayne  G.  Lee,  CE,  is  Commander  and  Director  of  WES.  Dr.  Robert  W. 
Whalin  is  Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC)  UNITS  OF  MEASUREMENT 
Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (metric) 


units  as  follows: 

Multidv 

Bv 

To  Obtain 

Btu  (International  Table) 
pound  (mass)  degree 

Fahrenheit 

per 

4,186.8 

joules  per  kilogram 

Kelvin 

Btu  (International  Table) 
inch  per  hour  .  square 
degree  Fahrenheit 

inch  . 

20.7688176 

watts  per  metre  Kelvin 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins* 

feet 

0.3048 

metres 

inches 

0.0254 

metres 

kips  (force)  per  inch 

1.213659 

kilonewtons  per  metre 

miles  per  hour  (U.  S.  statute) 

1.609347 

kilometres  per  hour 

pounds  (force)  per  square 

inch 

6.894757 

kilopascals 

pounds  (mass)  per  cubic  inch 

27,679.899 

kilograms  per  cubic  metre 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

*  To  obtain  Celsius  (C)  temperature  reading  from  Fahrenheit  (F)  readings,  use 
the  following  formula:  C  -  ( 5/9 )(F-32).  To  obtain  Kelvin  (K)  readings,  use 
K  -  ( 5/9 ) ( F -  32  )  +■  273.1  5. 


THERMAL  STRESS  ANALYSES  OF  MISSISSIPPI 
R.VER  LOCK  AND  DAM  26 (R) 

PART  I:  INTRODUCTION 
Background 

1.  In  November  of  1983,  the  Waterways  Experiment  Station  (WES)  was  asked 
by  St.  Louis  District,  Corps  of  Engineers,  to  analyze  the  lower  gate  monolith 
and  one  intermediate  monolith  of  Mississippi  River  Lock  and  Dam  26(R)  for 
thermal  and  construction  induced  stresses  and  also  normal  operating  conditions 
using  three-dimensional  (3-D)  finite  element  methods  (FEM) .  It  was  also  re¬ 
quested  that  WES  use  one  of  the  computer  progams  ADINA,  ANSYS ,  or  ABAQUS  to 
perform  the  analysis.  The  general  objective  of  this  analysis  effort  was  to 
adapt  a  general-purpose  finite  element  program  with  state-of-the-art  numerical 
formulations,  large  element  libraries,  and  easily  implemented  user-defined 
material  model  to  the  problem  of  thermal  stress  analysis  of  mass  concrete 
structures.  In  a  cooperative  effort  the  Concrete  Technology  Division  (CTD)  of 
the  Structures  Laboratory  (SL)  and  the  Information  Technology  Laboratory  (ITL) 
prepared  and  submitted  a  cost  estimate  and  scope  of  work  to  conduct  the  analy¬ 
sis.  Personnel  of  the  St.  Louis  District  (LMSED),  CTD,  and  ITL,  met  at  WES  in 
December  1983  to  review  and  discuss  the  cost  estimate  and  proposal.  All  pre¬ 
sent  at  this  meeting  agreed  to  the  general  analysis  approach  to  be  used,  but  it 
was  pointed  out  that  the  effective  solution  of  a  stress  analysis  problem  such 
as  this  might  require  significant  modifications  of  currently  available  general- 
purpose  finite  element  programs.  However,  it  was  anticipated  that  an  accept¬ 
able  solution  could  be  obtained  with  only  minor  modifications  of  existing 
finite  element  programs. 

2.  Monoliths  L-13  and  L-17  were  selected  for  the  analysis  so  that  two- 
dimensional  (2-D)  and  three-dimensional  (3-D)  effects  could  be  studied  both  in 


6 


terms  of  thermal  gradients  and  stress  distributions.  Symmetric  half -section 
isometric  projections  of  monoliths  L-13  and  L-17  are  shown  in  Figure  1. 


Ob i ective 

3.  The  objective  of  this  investigation  was  to  conduct  2-  and  3-D  ther¬ 
mal  stress  analyses  of  selected  monoliths  of  Lock  and  Dam  26(R)  using  a  modern 
general-purpose  finite  element  program. 


Scone 


4.  To  accomplish  the  objectives  of  this  study  within  the  guidelines  set 
by  the  St.  Louis  District,  a  four-phase  investigation  plan  was  developed.  The 
origionally  developed  plan  is  presented  below. 

Phase  _ Finite  Element  Investigation  Plan _ 

I  Selection  of  most  effective  finite  element  program(s)  to  accom¬ 

plish  the  goals  of  the  investigation.  This  involves  evaluating 
the  capabilities  of  the  available  three-dimensional  programs  to 
model  the  problem. 

II  Solution-method  verification  based  on  test  case  problems.  This 

involves  solving  simple  problems  for  which  solutions  are  avail¬ 
able  and  verifying  that  the  program  gives  exact  or  reasonable 
results  as  appropriate. 

Ill  Solution  of  two-dimensional  problems.  Two-dimensional  solutions 

will  be  obtained  for  monolith  L-13  using  the  WES2DT  program  which 
has  been  verified  based  on  previous  projects  where  2-D  approxi¬ 
mations  were  appropriate.  The  WES2DT  program  also  has  a  simple 
creep  model.  Two-dimensional  solutions  will  also  be  obtained 
for  monolith  L-13  using  the  new  code  and  compared  with  results 
from  the  WES2DT  program.  A  decision  was  made  later  to  include 
2-D  analyses  of  monolith  L-17  with  both  the  new  program  and 
WES2DT.  Results  from  these  2-D  analyses  form  a  basis  for 
interpreting  and  evaluating  results  from  3-D  analyses. 

IV  Solution  of  three-dimensional  problems.  A  3-D  analysis  will  be 

conducted  on  monolith  L-13.  This  analysis  is  primarily  designed 
to  demonstrate  the  effects  of  3-D  thermal  gradients  on  essen¬ 
tially  a  2-D  structural  geometry.  A  fully  3-D  analysis  will  be 
conducted  on  monolith  L-17  which  is  truly  a  3-D  structure . 
Although  at  the  beginning  of  the  study  it  was  intended  that 
complete  3-D  thermal  stress  analyses  would  he  conducted,  the 
computer  costs  and  limitations  on  computer  resources  prohibit t  •! 


PART  II:  SELECTION  OF  THE  FINITE  ELEMENT  PROGRAM 


General 


5.  The  approach  taken  in  selecting  the  finite  element  program  to  be 
used  in  the  Lock  and  Dam  26 (R)  analysis  was  to  review  the  capabilities  and 
flexibility  of  each  recommended  program  in  light  of  the  requirements  of  the 
problem  to  be  analyzed.  Specific  requirements  included  capability  to  simulate 
incremental  construction,  a  large  element  library  with  capability  to  implement 
user-defined  material  models  with  relative  ease,  capability  to  model  signifi¬ 
cant  numbers  of  reinforcing  bars  with  relative  ease,  efficient  numerical  solu¬ 
tion  procedures  with  flexible  means  for  selecting  solution  time  steps,  etc.  In 
addition,  it  was  desired  that  the  program  be  user-oriented,  receive  a  high 
caliber  of  technical  and  scientific  support  from  the  developer,  and  have  a  high 
potential  for  staying  at  the  state-of-the-art  level. 

6.  The  review  of  the  FEM  programs  was  conducted  by  personnel  of  the  ITL 

and  the  CTD.  The  review  consisted  of  discussing  actual  experiences  with  the 
different  programs,  reviewing  technical  journal  articles  which  compare  FEM 
programs  (1)  ,  meeting  with  representatives  of  Cybernet  Division  -  Control  Data 

Corporation,  and  discussions  with  personnel  of  other  governmental  laboratories 
directly  involved  with  FEM  applications. 

7.  Based  on  the  criteria  and  review  method  discussed  above,  the  review 
committee  unanimously  agreed  to  select  ABAQUS  for  the  Lock  and  Dam  26 (R)  pro¬ 


ject. 


ABACI'S  -  General  Descrintior 


8.  ABAQLS .  developed  bv  Hibbitt,  Karlsson.  and  Sorensen .  Inc . .  is  a 
general  -  purpose  structural  and  heat  transfer  ana! vs  is  program  The  theoret ii 
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formulation  is  based  on  the  finite  element  stiffness  method  with  some  hybrid 
formulations  included  as  necessary.  The  program  includes  both  user  and  automa¬ 
tic  control  of  solution  step  size.  Input  is  in  free  format,  key  worded,  and 
makes  use  of  set  definitions  for  easy  cross  reference.  A  broad  element  library 
is  included  in  ABAQUS .  and  any  combination  of  elements  can  be  used  in  the  same 
model.  A  wide  variety  of  constitutive  models  is  also  provided  in  ABAQUS,  and 
these  model  can  essentially  be  used  with  any  element  type.  User-defined  mate¬ 
rial  models  are  incorporated  with  relative  ease  through  the  UMAT  subroutine. 
Reinforcement  (rebar)  can  be  added  to  any  element.  Static  and  dynamic  response 
in  stress  analysis  can  be  conducted  as  well  as  steady-state  and  transient  heat 
transfer  problems.  The  incremental  construction  problem  can  be  effectively 
simulated  through  the  model  change  option  where  previously  defined  elements 
can  be  included  or  removed  from  the  analysis  in  a  specified  solution  step. 
Specific  details  of  element  type  and  boundary  conditions  used  will  be  presented 
in  the  various  analysis  phases  of  this  report. 
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PART  III:  ABAQUS  PROGRAM  -  CODE  VERIFICATION  AND  MODIFICATION 

General 

9.  The  code  verification  phase  was  to  serve  two  purposes:  (a)  to 
verify  the  operation  of  the  program  on  problems  with  known  solutions  and  its 
extension  to  problems  where  results  were  not  known  but  could  be  analyzed  by  the 
investigators  to  determine  that  solutions  followed  rational  and  expected 
behavior  patterns,  and  (b)  to  verify  the  investigators'  ability  to  direct  the 
program  operation  to  solve  desired  problems  by  using  the  many  different  program 
features . 


Code  Verification 

10.  To  verify  the  code,  several  simple  example  problems  were  used  and 
are  discussed  in  terms  of  two  groups.  Group  A  consists  of  problems  with  known 
solutions.  The  finite  element  (FE)  grids  for  these  example  problems  are  shown 
in  Figure  2.  Group  B  consists  of  the  example  problem  shown  in  Figure  3  which 
was  used  to  exercise  and  test  most  of  the  procedures  needed  to  solve  the  mass 
concrete  incremental  construction  problem  under  consideration. 

Group  A:  problems  with  known  solutions 

11.  Figure  2a  shows  a  square  plate  with  constant  temperature  specified 
on  each  boundary.  The  steady-state  temperature  distribution  was  determined 
with  ABAQUS  and  checked  against  the  closed  form  solution.  Very  good  agreement 
was  found  between  these  results. 

12.  Figures  2b  and  2c  show  two  s ingle - e 1 ement  models  using  the  eight- 
node  plane  strain  element.  These  models  were  subjected  to  two  temperature 
steps  of  10  degrees  F  each  while  E  changed  with  age  All  nodes  were 

res 'rained  during  the  temperature  change  in  model  2b.  whereas,  onlv  two  nodes 
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Figure  2. 


Finite  element  grids  used  in  ABAQUS  code  verification  for  Group  A 
problems  with  known  solutions. 
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were  partially  restrained  in  model  2c.  Results  showed  correct,  equal  stresses 
in  the  x  ,  y  ,  and  z  normal  directions  at  each  calculation  step  proportional 
to  the  temperature  change  and  E  values.  All  nodes  except  the  two  shown  were 
free  to  translate  in  model  2c.  The  results  showed  that  the  restraint  provided 
in  the  z  normal  direction  by  the  plane  strain  conditions  produced  correct, 
equal  stresses  in  the  x  and  y  directions  due  to  the  Poisson  effect  at  each 
calculation  step. 

13.  Figure  2d  shows  a  single-element  model  using  the  eight-node,  2-D 

heat  transfer  element.  In  this  example  problem,  the  user  subroutine  DFLUX  was 

tested.  DFLUX  was  wruten  to  provide  concrete  heat  of  hydration  input  to 

ABAQUS  thermal  simulations.  The  adiabatic  boundaries  of  the  model  prevented 

heat  exchange.  The  results  showed  that  the  sequences  of  temperature  changes  in 

the  model  exactly  reproduced  the  heat  generation  data  provided  to  DFLUX. 

Group  B:  problems  simulating 
incremental  construction 

14.  Figure  3a  shows  a  problem  consisting  of  24  eight-node  elements 
which  was  used  to  simulate  most  of  the  procedures  needed  in  the  actual  thermal 
stress  analyses  of  monoliths  L-13  and  L-17.  This  example  has  a  soil  foundation 
and  two  lifts  of  concrete  which  were  placed  at  3-day  intervals.  The  tempera¬ 
ture  distribution  is  determined  in  the  FE  system  throughout  the  incremental 
construction  process.  Next,  the  grid  in  Figure  3b  was  used  to  perform  the 
stress  analysis.  The  soil  foundation  was  not  included  in  the  stress  analysis, 
but  linear  springs  to  simulate  piles  were  added  to  support  the  structure  just 
as  in  the  actual  problem.  The  structure  was  then  constructed  in  two  lifts  to 
simulate  the  incremental  construction  process.  The  temperatures  which  were 
calculated  in  the  earlier  temperature  analysis  were  included  as  loads  to  the 
stress  analysis.  The  following  procedures  were  examined. 

a.  Eight-node,  2-D  element  behavior  for  temperature  analysis  and 
stress  analysis. 
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Figure  3. 


(b)  Grid  used  for  stress  analysis  during  incremental  con¬ 
struction  with  temperature  loadings 

I.  Finite  element  grids  used  in  ABAQUS  code  verification  for  Group  B 
problems  simulating  incremental  concrete  construction. 


b.  Linear  spring  element  behavior. 

c.  Gravity  loading. 

d.  Thermal  stresses. 

e.  Incremental  Construction  behavior. 

f.  Combination  of  temperature  analysis  and  stress  analysis. 

g.  Time -dependent  modulus  of  elasticity. 

h.  Gap  element  behavior. 

i.  Applied  pressures  and  concentrated  loads. 

j.  Plotting  of  displaced  grid  and  stress  contours. 

15.  During  the  analysis  of  the  results  of  this  example  it  was  deter¬ 
mined  that  the  incremental  construction  process  was  not  simulated  properly. 

The  program  developers  were  asked  to  supply  a  version  of  ABAQUS  which  could 
perform  this  task  as  required.  To  obtain  this  capability  the  new  ABAQUS 
Version  4.5  was  developed. 

16.  To  model  the  initial  temperature  equilibrium  behavior  at  the  inter 
face  of  an  intact  element  and  a  newly  placed  element,  the  use  of  gap  elements 
was  examined.  After  making  runs  without  elements,  with  gap  elements,  and  with 
modified  temperatures  at  those  interfaces,  it  was  determined  that  adequate 
results  were  obtainable  by  using  modified  temperatures  along  only  the  interfac 
between  the  soil  foundation  and  the  first  lift  of  concrete  during  the  tempera¬ 
ture  analysis.  This  was  very  desirable  since  the  use  of  gap  elements  would 
have  greatly  complicated  that  grid  generation  for  the  actual  models  to  be 
analyzed . 


17.  As  a  result  of  the  verification  examples,  four  modifications  were 
made  to  the  program.  The  first  modification  was  to  enable  the  program  to  use 
time -dependent  material  properties  (  E  and  v  ) .  The  second  modification  was 
to  enable  the  program  to  perform  the  incremental  construction  analysis  in  a 
more  consistent  manner.  Although  the  program  was  already  able  to  remove  and 
include  elements  of  the  structure  as  desired,  this  procedure  had  to  be  modified 
to  be  consistent  with  the  E  and  v  time  variations  for  newly  added  mate¬ 
rials.  This  was  implemented  in  ABAQUS .  The  third  required  modification  to 
ABAQUS  was  to  develop  and  implement  a  time  -  dependent ,  aging  creep  model  for 
material  behavior.  This  was  performed  by  adding  a  new  material  subroutine  to 
the  program.  This  subroutine  has  been  developed  and  implemented  in  ABAQUS  for 
2-D  problems,  and  a  version  for  3-D  problems  is  being  developed.  The  fourth 
modification  was  to  produce  the  user  subroutine  DFLUX  to  permit  incremental 
concrete  heat  of  hydration  data  to  be  generated  by  ABAQUS  for  the  incremental 
thermal  calculations. 


PART  IV:  WES2DT  PROGRAMS 


General 

18.  This  section  describes  the  two-dimensional  WES2DT  FE  programs  that 
have  been  used  at  WES  for  more  that  10  years  for  conducting  thermal  studies  of 
Corps  projects.  WES2DT  is  based  on  two  FE  thermal  analyses  programs  prepared 
for  the  Corps  during  the  late  1960's.  The  first  program,  developed  by  Dr. 
Edward  Wilson  of  the  University  of  California  at  Berkeley  (2)  and  modified  for 
use  at  the  WES,  calculates  the  temperatures  within  a  mass  concrete  structure. 

A  second  program,  written  by  R.  S.  Sandhu  et.  al .  also  at  Berkeley  (3)  and 
modified  at  WES,  calculates  the  thermal  stresses  and  strains  within  the  struc¬ 
ture  resulting  from  gravity  and  the  thermal  loads.  In  the  WES2DT  program 
construction  layer  (lift)  and  material  interfaces  must  correspond  to  an  element 
boundary.  The  element  type  used  in  WES2DT  is  a  four-node  quadrilateral  with 
linear  temperature  and  displacement  fields. 

Temperature  Calculation  Program 

19.  The  temperature  program  calculates  temperatures  at  each  node  in  the 
FEM  model.  These  calculations  are  based  upon  concrete  placement  temperature, 
heat  generated,  and  the  thermal  properties  of  the  concrete  which  govern  heat 
flow  within  an  element  and  loss  or  gain  across  boundaries  due  to  ambient 
conditions  which  are  controlled  by  a  surface  heat  transfer  coefficient  simula¬ 
ting  wind  or  surface  insulation  or  both.  Calculated  temperatures  are  output  at 
prescribed  intervals  for  all  nodes  in  the  model  at  a  particular  stage  of 
construction.  A  value  of  temperature  is  determined  for  each  new  concrete 
element  at  6  hr  after  placement  which  is  assumed  to  be  the  final  temperature  at 
which  an  element  is  stress  free.  Element  stress- free  temperatures  and  calcu¬ 
lated  nodal  temperatures  are  then  used  for  stress /s tra in  calculation  by  WF.S2DT . 
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Stress  and  Strain  Calculation  Proeram 


20.  This  program  calculates  the  displacements  at  each  node  and  the 
strains  and  stresses  developed  in  each  element  in  the  FEM  model  due  to  thermal 
and  gravity  loads.  It  was  discovered  that  the  program  was  designed  with  a 
plane  strain  formulation  in  which  both  z-direction  stresses  and  strains  are 
zero.  In  thermal  stress  problems  this  eliminates  z-direction  Poisson  effects 
on  x-  and  y-direction  stresses.  When  creep  is  considered,  stresses  at  each 
time  step  in  the  analysis  are  modified  for  stress  relaxation  at  constant 
strain.  Creep  parameters  are  stored  and  the  change  in  stress  stored  as 
residual  stress  to  be  included  in  the  next  time  step  analysis.  When  these 
stored  values  are  applied  during  the  next  time  step  analysis,  strains  are  then 
modified  for  creep. 

21.  A  modification  to  account  for  pile  restraint  was  incorporated  into 
the  thermal  stress  program  as  a  direct  result  of  requirements  for  previous 
thermal  studies.  The  pile  element  used  is  a  simple  two-node  element  modeled 
mathematically  as  a  linear  spring.  Actual  pile  stiffness  data  are  used  to 
determine  stiffness  of  the  pile  elements.  The  individual  stiffnesses  in  a  pile 
group  are  totaled  so  that  the  stiffnesses  can  be  applied  on  average  horizontal 
area  bases  rather  than  a  discrete  pile  basis  to  facilitate  2-D  analysis  with 
piles  included.  The  pile  groups  in  this  study  were  rows  of  piles  parallel  to 
the  flow  axis  of  each  monolith. 

22.  A  one  -  dimens ional  bar  (truss)  element  is  used  to  simulate  rein¬ 
forcement  steel.  Since  the  program  is  2-D,  the  reinforcement  parallel  to  the 
model  is  input  as  the  equivalent  cross-sectional  area  per  unit  depth.  Rein¬ 
forcement  and  pile  elements  are  superimposed  on  the  finite  element  model 
sharing  nodes  with  quadrilateral,  plain  strain  elements. 
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PART  V:  INPUT  DATA  AND  MODELED  PARAMETERS 


General 

23.  This  section  describes  the  input  data  and  modeled  parameters  used 
in  all  FEM  analyses  conducted  for  this  investigation.  The  first  section 
describes  the  materials  properties  for  concrete  and  foundation  as  well  as  pile 
stiffnesses  used.  The  second  section  describes  the  construction  parameters  and 
boundary  conditions  assumed  in  the  analyses. 


Concrete  and  Other  Properties 


24.  The  input  data  used  for  the  thermal  and  mechanical  properties  of 
concrete  were  either  the  same  values  used  in  earlier  thermal  studies  for  Lock 
and  Dam  26(R)  or  were  based  upon  tests  conducted  in  the  interim  period  since 
these  studies  were  completed.  The  concrete  is  assumed  to  have  a  nominal  com- 
pressive  strength  of  3000  psi  using  Type  II  cement  with  a  heat  of  hydration 
limit  and  25  percent  replacement  by  solid  volume  of  pozzolan  (fly  ash). 

25.  Concrete  properties  data  carried  over  from  previous  studies  in¬ 
cluded  adiabatic  temperature  rise  and  creep.  New  data  based  upon  tests  of  Lock 
and  Dam  26(R),  Phase  I,  project  mixture  4A  includes  modulus  of  elasticity, 
Poisson's  ratio,  specific  heat,  thermal  conductivity,  and  coefficient  of  linear 
thermal  expansion. 

Age -depr  ndent  concrete  properties 

26.  Adiabatic  temperature  rise  data  assumed  in  the  study  are  shown  in 

Figure  4a  Modulus  of  elasticity  as  a  function  of  age  used  bv  the  WES  FEM 

program  WF.S2DT  and  the  ABAQUS  UMAT1  subroutine  (without  creep)  is  shown  in 
Figure  4b. 

*  A  table  of  factors  converting  non-SI  units  of  measurement  to  SI  (metric’) 
units  is  presented  on  page  5. 
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27.  Creep  properties  used  in  the  WES2DT  code  were  taken  from  tests  conducted  at 


WES  on  concrete  from  Port  Allen  Lock  (4).  Generally,  creep  strains  for  this 


concrete  are  low,  thus  conservative,  for  mixtures  in  the  current  investigation.  The 


creep  data  were  fit  to  McHenry's  equation  (5,6), 


cc(<r,t,T)  -  c r£  Ak(T)  (l-e-mi^'T) ) 
i-1 


where 


c r-  applied  stress 


creep  strain 


t  -  time  after  placement 


T  -  age  at  loading. 


N  -  2  was  found  to  give  a  satisfactory  fit  of  experimental  data.  Values  of 


creep  relaxation  coefficients  Al  and  A2  versus  time  are  given  in  Figure  4c. 


Values  of  constants  -  0.45  and  m2  -  0.0285  were  used.  Figure  4d  shows  the 


creep  data  at  several  loading  ages  produced  by  the  creep  coefficients  shown  in 


Figure  4c  using  McHenry's  equation. 


28.  Creep  modeling  in  ABAQUS  was  accomplished  with  an  aging  creep  model 


that  was  developed  in  a  related  research  program  at  WES  (7).  The  2-D  aging 


creep  model  with  cracking  was  developed  in  the  ABAQUS -UMAT  subroutine  format. 


The  model,  designated  UMAT 2  for  this  investigation,  includes  the  effects  of 


Jgi-ng  on  the  elastic  modulus  and  cracking  strength  and  the  effects  of  changing 


temperatures  on  the  creep  compliance,  the  elastic  modulus,  and  the  ultimate 


cracking  strength.  Creep  properties  are  given  in  the  form 


cc/cT  -  [A(l-e'rC)  -t-  Bt  1  e' 


where 


e  -  creep  strain  due  to  stress  o 


3"  -  applied  stress 
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FUNCTION  OF  AGE  AS  USED  IN  THE  THERMAL  STRESS-STRAIN  RUNS 


LOADING  AGES 


IGiJRE  Ad.  SPECIFIC  CREEP  OF  CONCRETE  PS  A  FUNCTION  OF  AGE  AT  LOADING 
AS  USED  IN  THE  THERMAL  STRESS-STRAIN  RUNS.  DATA  COMPUTED 
BY  MCHENRY ' S  EQUATION  USING  COEFFICIENTS  FROM  FIGURE 
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T  -  temperature 
R  -  gas  constant 
r,A,B,Q  -  material  constants. 

The  elastic  modulus  and  ultimate  strength  can  be  expressed  as  functions  of  age 
(t)  and  temperature  (T)  as  E(t,T)  and  ou(t,T),  respectively.  The  cracking 
strain  (  cc  )  is  assumed  to  be  10  percent  of  the  absolute  value  of  the  compres 
sive  strain  at  ultimate  strength.  Poisson's  ratio  (  v  )  is  assumed  to  be 
constant.  The  creep  properties  and  the  properties  E(t,T),  ou(t,),  e^,  and  v 
are  included  in  a  separate  subroutine  in  a  form  that  can  be  modified  by  the 
user.  The  E  versus  time  data  used  in  UMAT2  was  based  upon  the  data  shown  in 
Figure  4b  except  that  the  numerical  fit  provided  by  the  model  developer  pro¬ 
duced  higher  values  of  E  for  the  first  two  days  of  age.  The  creep  data  used 
came  from  tests  of  silica  fume  concrete  and  produced  specific  creep  values 
slightly  higher  than  the  creep  data  shown  in  Figure  4d.  Shrinkage  data  was 
linear  until  1  day  at  a  rate  of  100  millionths  per  day  and  non-linear  there¬ 
after  with  typical  values  of  approximately  300,  350,  and  400  millionths  at  10. 
25,  and  50  days,  respectively. 

Other  properties 

29.  Additional  properties  data  used  as  input  to  the  computer  programs 
are  shown  in  Table  1.  Included  are  data  for  concrete,  the  soil  foundation,  an< 
air  when  modeling  the  internal  voids  of  monolith  L-17.  The  values  for  some  of 
the  properties  are  modified  to  accomplish  a  modeling  technique.  For  example, 
the  densit”  and  modulus,  E,  of  the  soil  in  the  VES2DT  sttess  analvses  were 
assigned  very  low  values  to  effectively  eliminate  the  soil  t ecause  "he  actual 
elements  could  not  be  removed.  An  abnormallv  high  value  of  thermal  conduct¬ 
ivity  was  assigned  to  air  to  simulate  the  combined  heat  transfer  effects  of 
t ion  and  conduction  when  air  elements  are  actual Iv  included  in  the 
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Table  1. 


Materials  Data  Used  In  FE  Investigation . 


I 

A 

j! 

v 


Concrete 


Thermal  conductivity  0.09789  Btu- in. /hr- in.  -°F 

Specific  heat  0.21  Btu/lb-°F  „ 

Density  0.08714  lb/in.  (151  lb/ft 

Coefficient  of  thermal  expansion  4.5  x  10*&  in./°F 
Poisson's  ratio  0.17 


s 


Foundation 

Thermal  conductivity  0.066  Btu-in./hr-in. ^-°F 

Specific  heat  0.45  Btu/lb-°F 

Density  (temperature  calculation)  0.0758  lb/in.  (131  lb/ft3) 
Density  (stress  calculation)  0.0000006  lb/in.  (or  -  0)* 

Coefficient  of  thermal  expansion  4.5  x  10*^  in./°F* 

Poisson's  ratio  0.35 

Modulus  of  elasticity  1.0  psi* 

Air  (for  voids  in  L-17  with  ABAQUS  temperature  calculations  only) 


Thermal  conductivity 
Specific  heat 
Density 


1000.  Btu-in./hr-in. ^-°F** 
0.24  Btu/lb-°F 
0.000046  lb/in.3 


*  Values  used  in  runs  with  WES2DT  program  only  to  eliminate  soil 
effects.  Foundation  not  included  in  stress  runs  with  ABAQUS. 

**  Modeled  value  for  air  used  combines  the  effects  of  convection  and 
conduction  in  an  enclosed  void  with  ABAQUS. 


Pile  Stiffnesses 


30.  Individual  pile  stiffnesses  for  H-piles  were  supplied  by  St.  Louis 
District.  Pile  stiffnesses  for  strong  soil  pile  support  were  used  because 
earlier  thermal  studies  on  several  projects  conducted  at  WES  indicated  that 
the  strong  pile  constant  assumptions  produced  the  highest  thermal  stresses. 

31.  Each  two-dimensional  FE  model  of  monoliths  L-13  and  L-17  represents 
a  plane  transverse  to  the  flow  axis  of  the  structure.  The  pile  stiffnesses 
actually  used  represent  the  average  stiffness  of  each  row  of  piles  parallel  to 
the  flow  axis  per  inch  thickness  into  the  model.  Figures  5a  and  5b  show  the 
plan  pile  layouts  for  monoliths  L-13  and  L-17. 
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32.  Monolith  L-13  In  monolith  L-13  there  are  either  8  or  15  piles  in 
the  flow  dimension  of  the  monolith  at  increasing  horizontal  distances  from  the 


axis.  All  individual  piles  were  specified  to  have  the  same  pile  stiffnesses 
and  all  are  oriented  with  the  strong  axis  in  the  horizontal  direction  of  the 
model  (transverse  to  the  flow).  Therefore,  the  stiffnesses  used  in  the  2-D 
analyses  differ  by  the  two  pile  densities  only.  These  stiffnesses  and  bases 
for  calculation  for  monolith  L-13  are  shown  in  Tables  2  and  3. 


Table  2. 

Individual  Pile  Stiffnesses.  Monolith  L-13 
(Strong  Pile  Constants) 


Fx  (horizontal, 
F  (horizontal, 
tz  (axial) 


strong  axis) 
weak  axis) 


Stiffness/Pile 
Ckips/in, ) 
54.6 
44.2 
896.4 


Table  3. 

Pile  Stiffnesses  Used  in  Two-Dimensional 
Stress  Analyses.  Monolith  L-13 
( Flow- direction  distance  *»  84  ft) 
(Stiffnesses  are  per  1.0-in.  thickness) 


Distance  of  Flow- Di rec t ion  Row  of  Piles 

_  from  Flow  Axis  (ft) _ 

0.  5.  10  15.  25.  35,  45.  55  65.  70.  75,  80 

No.  of  Piles  in  Row  15  8  15 

Stiffness  (k) 

(kips/in. ) 

Vertical  13.349*  7.1143  13.339 

Horizontal  0.8125*  0.4333  0.8125 


*  Values  at  0  ft  reduced  by  (  0.5  x  k  )  due  to  location  at  symmetric  boundary 


33.  Monolith  L-17  In  monolith  L-17  there  are  either  11  or  22  piles  in 
the  flow  direction.  The  transverse  distance  from  the  flow  axis  outward  was 
divided  by  St.  Louis  District  into  three  pile  stiffness  zones  as  shown  in 
Figure  5b.  Piles  in  zones  one  and  two  have  the  same  individual  pile 
stiffnesses,  however,  the  strong  axis  of  the  piles  in  zone  one  are  transverse 
to  the  flow  axis  while  the  strong  axis  of  piles  in  zone  two  are  in  the  flow 
direction.  Zone-three  piles  are  also  oriented  with  the  strong  axis  in  the  flow 
direction.  The  stiffnesses  used  and  bases  for  calculation  for  monolith  L-17 
are  shown  in  tables  4  and  5.  The  stiffnesses  were  determined  for  each  flow- 
direction  row  of  piles  individually  and  not  for  an  area  average  of  all  piles  in 
a  pile  zone. 

34.  In  3-D  models  of  monoliths  L-13  and  L-17,  individual  piles  are 
discretely  modeled.  Therefore,  the  stiffness  values  used  are  those  provided  in 
the  first  table  above  for  each  monolith.  Piles  were  oriented  in  the  directions 
described  earlier  and  shown  in  Figure  5. 


Table  4. 

Individual  Pile  Stiffnesses.  Monolith  L-17 
(Strong  Pile  Constants) 


Stiffness/Pile 
(kiDs/in . ) 

Zones  1  and  2 
(HP  14  x  73) 

Zone  3 

''HP  14  x  11") 

(horizontal , 

strong  axis) 

66 . 88 

66 . 88 

(horizontal , 

weak  axis) 

44 . 6 1 

36  .  10 

(axial) 

2375. 

►— * 
t' 

00 

* 


As  per  axis  direction  convention  provided  with  data  fr''m  St.  Louis 
District . 


Table  5 . 


Pile  Stiffnesses  Used  in  Two-Dimensional 
Stress  Analyses.  Monolith  L-17 
(Flow-Direction  distance  -  116  ft) 
(Stiffnesses  are  per  1.0-in.  thickness) 


Distance  of  Flow-Direction  Row  of  Piles 
_ from  Flow  Axis  (ft) _ 


15, 

20, 

25, 

70, 

75, 

30, 

35, 

40, 

60, 

80, 

85, 

0 

5 

10 

45. 

50__ 

55 

65 

90 

Pile  Zone 

3 

3 

3 

2 

1 

1 

No.  of  Piles/Row 

22 

22 

11 

11 

11 

22 

Strong  Axis 

with 

with 

with 

with 

with 

transverse 

Direction 

flow 

flow 

flow 

flow 

flow 

to  flow 

Stiffness  (k) 
(kios/in. ) 

Vertical 

11.680* 

23.359 

11.680 

18.768 

18.768 

37 .536** 

Horizontal 

0.285* 

0.570 

0.285 

0.352 

0.528 

1.057** 

*  Values  at  0  ft  reduced  by  (  0.5  x  k  )  due  to  location  at  symmetry  boundary. 

**  Values  at  90  ft  reduced  by  (  0.663  x  k  )  due  to  larger  area  of  influence 
associated  with  this  row  of  piles. 


Construction  Parameters  and  Boundary  Conditions 


Lift  heights 

35.  Lift  heights  used  in  the  FF.  models  for  both  monoliths  L-  1  3  and  L-17 
were  based  upon  those  provided  in  the  plans  supplied  hv  the  District.  Five-ft 
lifts  were  used  in  the  base  and  8-  to  10- ft  lifts  were  used  in  the  wall  of 
monolith  L-13.  Lift  heights  ranged  from  3  ft  to  5  ft  in  the  base,  6  ft  in  The 
lower  wall,  and  5  to  5.5  ft  in  the  upper  wall  of  monolith  L-l'.  The  actual 
locations  of  lifts  are  shown  on  the  FF  models  in  Figures  7a-  'e  in  Part  VI  of 
this  report . 
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Placement  rates 

36.  A  concrete  placement  interval  of  5  days  between  lifts  was  used 

throughout  the  study.  Although  actual  placement  may  progress  more  slowly,  a 

rate  of  5  days  per  lift  was  used  as  the  most  possible  case.  This  rate  also 

allows  less  time  for  lift  surfaces  to  cool  between  placements  of  lifts  than  at 

slower  rates  of  placement.  The  higher  internal  temperatures  and  thermal 

stresses  that  result  become  the  maximum  values  attainable  at  the  most  rapid 

allowable  placement  rate. 

Ambient  air  temperature 
and  foundation  temperature 

37.  Air  temperatures  used  were  provided  by  the  District  and  are 
representative  of  the  project  site.  Figure  6  contains  the  mean  ambient 
temperature  versus  time  of  the  year  data  used  in  the  study.  The  data  are 
characterized  by  a  sinusoidal  curve  peaking  in  August  and  having  minimum 
temperatures  in  January.  Including  daily  variation  in  temperature  was  felt  to 
be  beyond  the  scope  of  this  study. 


38.  Foundation  temperature  was  assumed  to  be  a  constant  55°  F  at  a  depth 
of  20  ft  below  the  lowest  elevation  of  each  monolith.  In  this  manner  the 
constant  temperature  was  specified  at  el.  344  ft  in  monolith  L-13  and  el.  338  ft 
in  monolith  L-17.  The  value  of  55°  F  represents  mean  annual  temperature  for  the 
project  site.  The  initial  temperature  distributions  within  the  soil  were 
determined  in  a  preliminary  finite  element  calculation  in  which  the  upper 
surface  of  the  soil  was  exposed  to  the  ambient  temperature  data  shown  in  Figure 
6  and  the  constant  value  of  55°  F  along  the  lower  model  bound,  ry.  This  heat 
flow  simulation  was  run  over  a  scaled  time  period  exceeding  one  year  so  that 
final  values  accounted  for  the  thermal  momentum  of  the  soil  mass  in  responding 
to  ambient  temperature  changes. 
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FIGURE  6  -  MEAN  DAILY  AMBIENT  TEMPERATURES  FOR  ALTON,  IL.  USED 

IN  TEMPERATURE  RUNS.  DATA  REPRESENTS  EXPECTED  NORMAL 
TEMPERATURES. 


Lacement  temperature 


39.  All  computer  runs  in  this  investigation  were  based  upon  a 


construction  start  date  of  1  July.  Previous  investigations  indicated  that  an 


early  summer  construction  start  generally  produces  highest  thermal  stresses 


when  compared  to  other  start  times  during  the  year. 


40.  A  maximum  placement  temperature  of  65°  F  was  used  in  all  areas  of 


the  structure.  During  summer  construction  this  will  probably  require  that  some 


of  the  mixing  water  be  added  as  crushed  ice. 


Thermal  boundary  conditions 


41.  The  lower  boundary  of  the  soil  was  fixed  at  55  F  as  described 


earlier.  No  horizontal  heat  flow  was  permitted  through  the  vertical  soil 


boundaries.  Heat  flow  was  also  not  permitted  horizontally  through  the  vertical 


centerline  of  the  monoliths  represented  by  the  extreme  left  side  of  the  models. 


42.  The  thermal  boundary  condition  that  controls  heat  exchange  between 


the  structure  and  the  ambient  air  is  a  surface  heat  transfer  coefficient.  It 


is  composed  of  a  convection  coefficient  that  defines  heat  exchange  with 


surrounding  air  as  well  as  a  conduction  heat  transfer  coefficient  which  defines 


the  heat  flow  through  the  formwork. 


43.  The  convection  heat  transfer  coefficient  is  based  upon  surface  air 


velocity.  An  air  velocity  of  10  mph  was  used  for  external  surfaces  and  a 


velocity  of  1  mph  was  used  for  the  surfaces  of  unenclosed  voids  such  as 


galleries,  culverts,  and  recesses.  It  was  assumed  that  unenclosed  voids  would 


be  exposed  to  and  ventilated  by  ambient  air,  but  with  reduced  air  flow.  The 


enclosed  voids  (internal  rooms)  in  monolith  L-17  were  treated  similarly  until 


the  enclosing  concrete  of  lift  8  was  placed.  At  void  closure,  the  FE  grid  in 


the  void  area  was  activated  as  air  elements.  This  permitted  the  temperature  of 


the  void  to  varv  as  a  function  of  surrounding  concrete  temperatures  as  ii 
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reality.  This  modeling  scheme  was  considered  adequate  even  though  convection 
heat  exchange  between  the  void  and  concrete  surfaces  and  the  vertical 
stratification  of  the  air  could  not  be  included. 

44.  On  external,  vertical  surfaces  and  non-enclosed  recesses  the 
insulating  effect  equivalent  to  3/4- in.  plywood  forms  was  included  in  the  heat 
transfer  coefficient  for  the  first  2  days  after  placement  of  a  lift.  After 

the  second  day  in  order  to  simulate  form  removal,  the  heat  transfer  coefficient 
included  only  convection  effects  of  the  wind.  On  internal,  vertical  formed 
surfaces  and  the  horizontal  upper  void  surfaces  of  monolith  L-13,  the  effect  of 
forms  was  left  in  place  until  seven  days  after  placement  of  the  enclosing  cover 
lift  to  simulate  the  forming  practices  in  such  areas.  In  monolith  L-17  all 
formwork  heat  transfer  effects  were  removed  at  the  time  that  the  void  cover 
lift  was  placed  as  described  in  Paragraph  43  above.  All  exposed  horizontal 
surfaces  were  given  heat  transfer  coefficients  equivalent  to  convection  wind 
effects  only.  No  supplemental  insulation  was  included.  In  3-D  runs,  the  flow- 
direction  end  of  the  monoliths  were  assumed  free.  That  is,  adjoining  monoliths 
were  not  assumed  to  exist  in  these  simulations. 

Mechanical  boundary  conditions 

45.  Because  it  was  assumed  that  the  piles  elements  would  carry  the 
concrete  loads,  the  soil  was  ignored  in  all  stress  runs.  The  symmetric 
boundaries  were  fixed  in  the  horizontal  directions.  All  other  surfaces  were 
assumed  free  including  the  ends  of  the  monolith  in  3-D  runs. 
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PART  VI:  FINITE  ELEMENT  MODELS 


46.  The  primary  goal  of  this  investigation  has  been  to  analyze  selected 
monoliths  of  Lock  and  Dam  26(R)  for  thermal  and  construction  induced  stresses 
based  upon  calculations  using  state-of-the-art  FE  technology.  While  efforts 
were  directed  toward  this  goal,  a  substantial  amount  of  developmental  work  and 
investigative  analyses  were  conducted  to  implement  the  mechanics  of  using 
ABAQUS  for  incremental  construction  thermal  stress  analysis,  to  assess  the 
requirements  for  acceptable  modeling  accuracy,  and  to  evaluate  the  results 

of  these  analyses.  These  factors  required  development  of  many  FE  models,  both 
two-  and  three-dimensional,  in  order  to  complete  this  investigation.  The 
following  sections  describe  these  models,  including  the  analytical  bases,  and 
the  other  considerations  taken  in  their  development. 

47.  Monoliths  L-13  and  L-17  were  requested  for  FE  analyses  in  this 
investigation.  L-13  represents  a  typical  lock  chamber  monolith,  which  due  to 
its  common  transverse  -  section  geometry,  was  considered  ideal  for  comparative 
two-  (2-D)  and  three-dimensional  (3-D)  analyses.  L-17  is  the  lower  gate 
monolith  which  is  more  massive,  and  has  geometric  features  which  include  large, 
non- symmetrically  located  internal  voids  or  rooms,  non- symmetrical  external 
recesses,  and  different  chamber  floor  thicknesses.  These  features  point 
strongly  toward  a  requirement  for  3-D  analysis.  Symmetric,  half-section 
isometric  projections  of  monoliths  L-13  and  L-17  are  shown  in  Figure  1. 
Internal  voids  in  L-17  are  indicated, 

48.  All  2-D  models  developed  for  this  investigation  represent  sectional 
planes  transverse  to  the  flow  axis.  The  monoliths  were  assumed  .symmetric 
relative  to  the  flow  axis,  therefore ,  each  2-D  model  consists  of  nm.--h.il  t  oi 
the  transverse  cross  section  In  all  cases  the  axis  of  .,vm:ret  rv  lies  at  the 
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lefc  side  of  the  model.  The  3-D  models  consist  of  quarter  sections  of  the 
monoliths  with  symmetric  axes  located  at  the  flow  axis  and  halfway  between  the 
upstream  and  downstream  faces.  Actually,  the  3-D  model  of  L-17  is  not  exact 
geometrically  because  the  monolith  is  not  symmetric  in  the  flow  direction. 

This  subject  is  further  explained  below. 

49.  In  all  FE  models  the  soil  was  included  for  the  temperature 
analysis.  With  ABAQUS ,  the  soil  was  removed  from  the  stress  analysis  with 
spring  (pile)  elements  added  to  support  the  structure.  The  WES2DT  program  did 
not  permit  removing  the  soil  elements.  As  an  alternative,  the  mechanical 
properties  of  the  soil  were  changed  to  appropriate  values  (see  Table  1)  which 
effectively  removed  the  soil.  Pile  elements  were  added  to  support  the 
structure  in  a  manner  identical  to  that  when  using  ABAQUS. 


50.  Based  primarily  on  past  experience  with  the  WES2DT  program,  a  grid 
was  selected  for  use  with  ABAQUS  that  was  expected  to  have  an  adequate  number 
of  degrees  of  freedom  to  capture  the  essential  features  of  the  temperature  and 
stress  response  expected  for  the  problem  especially  considering  the  higher 
order  elements  available  in  ABAQUS.  This  FE  grid  referred  to  as  the  coarse 
grid  is  shown  in  Figure  7a  and  is  characterized  by  two  elements  per  lift, 
vertically,  and  two  elements  wide  across  individual  sections  of  the  culvert  and 
wall  stem.  To  insure  that  this  grid  was  adequate,  a  more  refined  grid,  Figure 
7b,  was  also  developed  to  be  used  for  the  2-D  analyses.  The  refined  grid  is 
characterized  by  three  elements  per  lift,  vertically,  and  three  elements  wide 
across  individual  sections  of  the  culvert  and  wall  stem.  The  assumption  made 
was  as  follows:  if  the  results  of  analyses  using  the  two  grids  are  in  close 

agt  eement ,  the  amount  of  refinement  in  the  first  r.r  id  would  be  sufficient  since 


similar  grids  have  been  used  in  the  past  with  WES2DT.  This  being  the  case, 
preparation  of  only  one  grid  each  would  be  required  for  2-D  analyses  of  L-17 
and  for  all  3-D  analyses. 

WES2DT  model  selection 

51.  The  FE  grid  of  L-13  developed  for  use  with  the  WES2DT  program  was 
based  upon  the  coarse-grid  model  used  with  ABAQUS,  but  with  an  exception.  The 
WES2DT  program  uses  four-node  elements  with  linear  temperature  distribution  and 
linear  displacement  functions  while  ABAQUS  was  run  with  eight-node  elements 
using  quadratic  temperature  and  displacement  functions.  In  order  to  make 
reasonable  comparisons  between  the  two  programs  the  grid  developed  for  WES2DT 
was  designed  to  coincide  with  the  nodal  locations  in  the  ABAQUS  model.  This 
process  yielded  four  elements  in  the  WES2DT  grid  for  every  one  in  the  ABAQUS 
equivalent  grid.  The  FE  grid  used  to  model  L-13  for  the  WES2DT  program  is 
shown  in  Figure  7c. 


Two-Dimensional  Models.  L-17 


ABAQUS  model  selection 

52.  The  element  density  criteria  used  to  establish  the  2-D  grid  for  L- 
17  was  similar  to  that  used  for  the  coarse-grid  FE  model  of  L-13,  namely,  using 
two  elements  vertically  per  lift.  The  element  density  horizontally  generally 
conforms  to  the  locations  of  rows  of  piles  which  are  oriented  parallel  to  the 
flow  axis  and  also  to  the  location  of  internal  voids.  The  2-D  section  of  L-17 
chosen  was  that  located  at  Station  26  +  13  near  the  downstream  end  which 
contained  the  smallest  amount  of  voids.  This  section  was  expected  to 
experience  the  largest  temperature  rise  and  was  the  same  one  used  for  an 
earlier  temperature  study  using  the  WES2DT  program.  This  section  passes 
through  the  smaller  of  the  internal  voids  in  L-17  and  through  a  machinery 
recess  at  the  top  of  the  monolith. 
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0  -  indicates  pile  element  locations 
*  -  elevations  are  in  inches 


Figure  7c.  Monolith  13,  two-dimensional  finite  element 
grid  u:  ed  with  WES-2DT  programs. 


53.  Voids  and  recesses  were  not  gridded  in  all  previous  FE  models. 

With  ABAQUS ,  elements  can  be  activated  or  deactivated  dynamically  as  the 
incremental  construction  process  is  simulated.  This  feature  provided  an 
opportunity  for  a  new  approach  in  FE  grid  preparation  for  L-17.  A  sectional 
grid  which  represents  the  maximum  floor  and  wall  thicknesses  was  fashioned.  By 
locating  the  grid  to  coincide  with  recesses  and  voids  as  well  as  construction 
joints  between  lifts,  virtually  any  cross  section  of  L-17  could  be  modeled  by 
merely  deactivating  the  appropriate  elements  which  represented  voids  or  recesses. 
This  was  particularly  useful  when  dealing  with  the  internal  voids. 

54.  In  L-13  the  culvert  and  upper  access  gallery  pass  through  the 
entire  length  of  the  monolith.  It  was  assumed  when  modeling  the  thermal 
boundary  conditions  that  these  voids  would  be  ventilated  by  ambient  temperature 
air,  but  with  a  lower  velocity  than  that  on  exposed  external  surfaces.  In  L-17 
the  internal  voids  are  exposed  to  ambient  temperature  air  only  until  the 
enclosing  concrete  is  placed.  Except  for  doorways  to  the  adjacent  voids,  these 
voids  are  enclosed.  After  closure,  the  air  temperatures  in  the  voids  must 
become  a  functions  of  the  concrete  surface  temperatures  surrounding  each  of  the 
voids . 

55.  The  ABAQUS  2-D  model  of  L-17  shown  in  Figure  7d  contains  the  full 
sectional  grid  as  described  above  with  two  shaded  areas  included.  The  elements 
in  the  upper  shaded  area  represent  the  machinery  recess  and  were  deactivated 
permanently.  The  elements  in  the  lower  shaded  area,  which  represents  the 
internal,  closed  void,  were  deactivated  only  until  the  covering  lift  (No.  8) 
was  placed.  The  void  elements  were  then  activated  with  the  properties  of  air. 
Thereafter,  these  air  elements  were  an  integral  part  of  the  model.  This 
procedure  was  an  effective  means  of  modeling  the  thermal  boundary  of  the  void. 

The  only  limitation  to  this  modeling  technique  is  that  convection  heat  transfer 
will  in  reality  cause  air  in  the  void  to  stratify  with  the  warmer  air  rising  to 
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Figure  7d.  Monolith  17,  two-dimensional  finite  element  model,  ABAQl'S, 
program  (316  elements) 
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the  ceiling.  This  minor  error  is  acceptable  in  light  of  the  benefits  derived 
from  implementing  this  modeling  technique. 

WES2DT  model  selection 

56.  The  2-D  FE  grid  of  L-17  designed  for  use  with  WES2DT  program  for  an 
earlier  thermal  study  of  L  &  D  26(R)  is  shown  in  Figure  7e.  This  grid, 
therefore,  does  not  compare  exactly  with  the  one  prepared  for  use  with  ABAQUS 
as  was  the  case  for  L-13.  However,  if  the  four-element- for-one  technique  for 
designing  a  comparable  WES2DT  FE  grid  from  an  ABAQUS  grid  was  followed  for  L-17 
as  was  employed  for  L-13,  the  WES2DT  grid  generally  compares.  One  major 
difference  is  that  the  WES2DT  grid  only  extends  vertically  through  lift  13  (el. 

420.5  ft.)  whereas  the  ABAQUS  grid  extends  the  full  16-lift  height  of  the 
monolith  (el.  434.5  ft.).  Through  the  first  13  lifts,  the  WES2DT  grid 
consisted  of  over  900  elements  which  was  considered  a  practical  upper  limit. 

Because  the  primary  areas  of  interest  were  the  chamber  floor  and  lower  wall, 
the  top  lifts  were  not  gridded.  Instead,  the  top  surface  of  lift  13  was 
insulated  appropriately  at  the  scheduled  time  of  placing  lift  14  to  simulate 
the  correct  temperature  conditions  in  the  monolith.  Equivalent  gravity  loads 

for  lifts  14-16  were  simulated  by  applying  surface  pressure  to  the  top  of  lift  13  at 
the  appropriate  times  of  placement  of  the  subsequent  lifts. 

Three-Dimensional  Model.  L-13 

57.  As  described  earlier,  the  3-D  FE  model  prepared  for  L-13  is  a 
quarter  section  of  the  monolith.  This  model  shown  in  Figure  8a  is  made  up  of 
five  element  planes  in  the  flow  direction  each  with  a  sectional  grid  identical 
to  the  coarse-grid  2-D  model .  The  four  interfaces  planes  between  these  five 
element  planes  coincide  with  the  location  of  those  complete  rows  of  piles 
oriented  transverse  to  the  flow  direction.  Figure  5a  shows  a  plan  view  of  the 
quarter-section  L-13  3-D  model  superimposed  on  the  pile  lavout  with  vertical 
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element  boundary  projections  shown.  The  grid  was  designed  to  model  the  geometry 
of  the  monolith  with  a  minimum  number  of  elements  in  the  flow  direction  which 
were  adequate  to  provide  acceptable  results  and  have  element  boundaries 
coinciding  with  pile  locations.  The  20-node  3-D  elements  that  were  used  have 
nodes  located  at  the  element  corners  and  midway  between  corners.  As  seen  in 
Figure  5a,  most  of  the  piles  coincided  with  these  nodal  locations  in  the  model. 
For  those  that  did  not,  a  capability  within  ABAQUS  was  used  that  permited 
creation  of  extra  nodes  at  these  unassigned  pile  locations. 

Three-Dimensional  Model.  L-17 

58.  It  was  recognized  very  early  in  the  investigation  that  proper  3-D 
modeling  of  L-17  would  require  a  half-section  symmetric  representation  of  the 
monolith.  This  was  due  to  the  non- symmetric  geometry  in  the  flow  direction. 
However,  constraints  identified  by  the  developers  of  ABAQUS  and  those 
recognized  by  the  investigators  indicated  that  a  half-section  model  was  not 
feasible.  In  fact  it  was  recognized  that  even  a  minimally  discretized  quarter- 
section  model  would  severely  tax  available  computer  and  project  resources  due 
to  the  massiveness  of  the  monolith.  The  primary  constraint  was  a  practical 
upper  limit  of  approximately  1000  elements. 

59.  It  was  evident  that  several  trade-offs  would  be  necessary  for  the 
L-17  3-D  model.  The  construction  plan  for  L-17  indicated  concrete  placement  in 
16  lifts.  It  had  been  determined  from  the  ABAQUS  incremental  construction 
evaluations  and  L-13  2-D  grid  comparisons  that  at  least  two  elements  vertically 
per  lift  were  preferred  for  adequate  calculation  accuracy.  It  was  also  evident 
after  examination  of  the  geometric  features  that  a  minimum  of  six  elements  were 
required  in  the  flow  direction  in  order  to  model  the  geometric  features,  which 
included  voids  and  recesses,  and  to  include  the  actual  flow  direction  axis  of 
symmetry.  Based  upon  the  5-ft  hv  5-ft  spacing  of  the  basic  pile  location  grid, 


it  was  evident  that  the  maximum  size  of  a  3-D  element  in  the  floor  slab  was  10- 
ft  square  in  plan.  Smaller  elements  were  necessary  in  the  wall.  These  criteria 
were  used  to  lay  out  a  preliminary  grid.  The  resulting  grid  size  was  close  to 
2000  elements  or  nearly  double  the  practical  limit. 

60.  To  keep  within  the  1000-element  practical  limit  several  changes  had 
to  be  made.  First,  except  for  lift  5,  each  lift  was  gridded  with  only  one 
element  vertically  rather  than  two.  The  number  of  element  planes  in  the  flow 
direction  had  to  be  reduced  from  six  to  five.  In  order  to  properly  model  the 
monolith  and  discrete  piles,  it  was  necessary  for  element  boundaries  in  the  flow 
direction  to  correspond  to  every  other  transverse  row  of  piles.  Therefore,  the 
resulting  quarter- section  boundary  transverse  to  the  flow  is  not  located  at  the 
monolith  midpoint  in  the  flow  direction.  The  3-D  nodes  is  slightly  less  than 
one-half  the  monolith  dimension  in  the  flow  direction. 

61.  The  basic  3-D  grid  resembled  an  reversed  "L"  shape  for  the  floor  and 
wall  similar  to  that  used  for  the  2-D  L-17  model.  All  voids  and  recesses  were 
similarly  obtained  by  deactivating  the  appropriate  elements.  Figures  8b  and  8c 
are  isometric  projections  of  the  L-17  3-D  grid.  Figure  8b  shows  the  basic  grid 
with  external,  deactivated  elements  shaded  while  Figure  8c  shows  the  grid  with 
the  deactivated  elements  removed.  Figure  8d  shows  the  basic  sectional  element 
plane  and  che  five  elements  planes  with  the  internal  and  external  deactivated 


elements  shaded. 
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PART  VII:  THERMAL  ANALYSES  OF  MONOLITH  L-13 


Presentation  of  Result; 


Two-dimensional 
temperature  analysis 

62.  Temperature  analyses  of  monolith  L-13  were  performed  using  both  the 

ABAQUS  and  WES2DT  programs.  Figures  9a  thru  9d  show  results  from  the  ABAQUS 

program  using  FE  grid  1  (Figure  7a)  at  four  stages  of  construction.  Figures 

10a  thru  lOd  show  the  results  of  the  ABAQUS  program  using  FE  grid  2  (Figure  7b) 

at  four  stages  of  construction.  Results  from  the  WES2DT  program  using  a  FE 

grid  (Figure  7c)  functionally  equivalent  to  that  of  ABAQUS  FE  grid  1  are  shown 

in  Figures  11a  thru  lid  for  the  same  four  stages  of  construction  used  for 

presenting  ABAQUS  data. 

Three  dimensional 
temperature  analysis 

63.  The  ABAQUS  program  was  used  to  make  a  3-D  analysis  of  monolith  L-13 
using  the  FE  grid  shown  in  Figure  4.  These  results  are  shown  using  contour 
plots  (Figures  12a  -  12d  thru  15a  -  15d)  which  show  temperatures  on  faces  of 
elements  in  the  five  element  planes  (Figure  8)  of  the  structure  at  four  stages 
of  construction.  Four  of  the  element  planes  are  internal,  one  is  an  outer 
surface . 


Discussion  of  Results 

Two-dimensional 
temperature  analysis 

64.  The  ABAQUS  results  for  FE  grids  1  and  2  were  virtually  identical, 
it  was  evident  from  these  runs  that  the  refinement  in  grid  1  was  adequate  for 
the  temperature  analysis.  The  comparison  of  ABAQUS  results  to  WES2DT  results 


was  also  excellent.  The  only  minor  discrepancy  came  in  the  contour 
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afcer  lift  5  was  placed  (Figures  9c,  10c,  and  11c)  in  the  lift  5  concrete. 

Here  the  ABAQUS  results  show  a  larger  area  with  temperatures  above  95°  F  in 

lift  5.  This  seems  more  realistic  since  the  base  slab  exhibited  a  similar 

contour  pattern  at  the  the  same  time  after  placement  (Figures  9a,  10a,  and 

11a).  Only  localized  differences  were  seen  between  these  different  analyses, 

with  maximum  values  and  the  general  temperature  distribution  in  very  good 

agreement  at  all  stages  of  the  construction  sequence. 

Three-dimensional 
temperature  analysis 

65.  Upon  examining  the  temperature  contour  plots  of  the  base  slab, 
Figures  12b  thru  12c,  it  is  evident  that  the  temperature  distribution  was 
constant  along  the  monolith  in  the  direction  of  flow  (z-direction) .  This  was 
expected  since  the  problem  closely  approximates  2-D  assumptions  for  constant 
boundary  conditions  on  the  outer  surface  and  centerline.  Only  in  the  vicinity 
of  five  feet  of  the  outer  surface  did  any  variation  take  place.  This  is  shown 
in  Figure  12a.  Similar  response  is  shown  by  the  contours  for  all  other  stages 
of  construction  (Figures  13  thru  15).  The  contours  within  the  monolith  (away 
from  the  exposed  outer  surface)  also  match,  almost  identically,  the  ABAQUS  2-D 
plots  from  Figures  9  thru  10.  These  contours  show  that  monolith  L - 1 3 ,  for  these 
boundary  conditions,  is  a  2-D  problem. 
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Figure  9a.  Temperature  distribution  in  structure  5  days  after  lift  2 
is  placed  (from  analysis  using  first  grid) 
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Figure  10b .  Temperature  distribution  in  structure  5  days  after  lift  4 
is  placed  (from  analysis  using  second  grid) 
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Figure  lOd.  Temperature  distribution  in  structure  7  days  af 
is  placed  (from  analysis  using  second  grid) 
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igure  lib.  Temperature  distribution  in  structure  5  days  after  lift  4  is  placed  from  WES-2DT 
program. 
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Figure  12b.  Temperature  distribution  in  node  plane  5  of  structure  5  days 
after  lift  2  is  placed. 
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Figure  12c.  Temperature  distribution  in  node  plane  4  of  structure  5  days 
after  lift  2  is  placed. 


l.o.  VAUiC 
I  *«.50 1*01 
i  *7. 00C -01 
3  *7.  50C *01 
M  -•  00E»0l 

s  *0.S0C’0l 
s  *9  OOC  *01 

7  «9.  SOC *0 1 

a  *i  ooc -02 

9  *1  0SE*02 

10  •!  I OE  *02 

11  « I  I  SC  *02 


LIFT2  -  M 1 3  -  30  CRIO  MODEL  -  NODE  PLANE  2 

STEP  <t  INCPENENT  «  osaaus  VERSION  n-s-ivt 


Figure  12e.  Temperature  distribution  in  node  plane  2  of  structure  5  days 
after  lift  2  is  placed. 
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Figure  12f.  Temperature  distribution  in  the  centerline  of  structure  5  days 
after  lift  2  is  placed. 
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Figure  14a.  Temperature  distribution  in  outer  surface  of  structure  5  days 
after  lift  5  is  placed. 
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Figure  15d.  Temperature  distribution  in  node  plane  3  of  structure  7  days 
after  lift  9  is  placed. 


i.o.  vxiue 

1  •S.OOE^OI 

2  ♦7.00E*0l 

3  .750E.0I 
1  •B.OOE.OI 

s  <a.riOE'Oi 

6  •9.00E.0I 

7  *9. 50E*0I 
6  *  1 . 00E.02 
9  *  1 . 05E-*02 

10  ♦!  10E«02 

11  - I . ISE«02 


LIFT9  -  M 1 3  -  30  GRID  MQOEL 


STEP  18  INCftEMENT  10 


NODE  PLANE  2 


aanous  version  u -s- i u7 


Figure  15e.  Temperature  distribution  in  node  plane  2  of  structure  7  days 
after  lift  9  is  placed. 
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PART  VIII:  STRESS  ANALYSES  OF  MONOLITH  L-13 


66.  The  ABAQUS  program  was  used  to  perform  the  2-D  stress  analysis  of 

monolith  L-13  for  both  instantaneous  gravity  turn-on  of  the  entire  structure 

and  an  incremental  construction  build-up  sequence  of  the  structure.  Figures 

16a  and  16b  show  the  displaced  structure  and  maximum  principal  (tensile)  stress 

contours,  respectively,  for  the  gravity  turn-on  analysis  for  a  Young's  Modulus 

(E)  of  3.12  x  10^  psi.  Figures  17a  and  17b  give  the  same  results  for  E  -  4.80 

x  10^  psi.  These  modulus  values  correspond  to  the  ACI  modified  and  normal 

compressive  strength  of  the  concrete  at  28  days,  respectively.  Figures  18a 

thru  18d  and  19a  thru  19b  show  the  displaced  structure  and  maximum  principal 

stress  contours,  respectively,  for  four  four  stages  of  construction  using 

incremental  construction  sequencing.  This  incremental  construction  analysis 

uses  the  material  modulus  calculation  routine  UMAT1.  WES2DT  was  used  to  make 

the  same  incremental  construction  analysis.  Figures  20a  thru  20d  and  21a  thru 

21d  show  the  resultant  displacement  vector  plots  and  principal  stress  vector 

plots,  respectively,  from  this  analysis. 

Two-dimensional 

gravity  and  thermal  loading 

67.  The  ABAQUS  program  was  used  to  perform  the  2-D  stress  analysis 
including  gravity  and  temperature  loading  for  FE  grids  1  and  2  using  the 
modulus  routine  UMAT1.  The  displaced  structure  plots  from  grids  1  and  2  are 
given  in  Figures  22a  thru  22d  and  24a  thru  24d,  respectively.  The  maximum 
principal  stress  contours  for  both  grids  are  given  in  Figures  23a  thru  23d  and 
25a  thru  25d,  respectively.  The  new  modulus  subroutine  UMAT2  (concrete  aging 
creep  model)  was  used  with  ABAQUS  with  grid  1  to  make  the  analysis  when  E 
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did  not  depend  on  temperature  and  E  did  depend  on  temperature.  The  displaced 
structure  and  maximum  principal  stress  contours  are  given  in  Figures  26a  thru 
26d  and  27a  thru  27d  for  E  not  a  function  of  temperature.  Figures  28a  thru 

28d  and  29a  thru  29d  show  the  displacement  and  maximum  principal  stress  con¬ 

tours  for  E  as  a  function  of  temperature.  A  WES2DT  analysis  was  made  of  this 
same  gravity  and  temperature  loading.  The  resultant  displacement  vector  plots 
are  shown  in  Figures  30a  thru  30d.  The  principal  stress  vector  plots  are  shown 
in  Figures  31a  thru  31d. 

Two-dimensional  gravity  and 
thermal  loading  including  creep 

68.  ABAQUS  was  used  in  incremental  construction  analyses  of  the  struc¬ 

ture  using  FE  grid  1  with  the  modulus  subroutine  UMAT2 .  Displaced- structure 
plots  are  shown  in  Figures  32a  thru  32d.  Maximum  principal  stress  contours  are 
shown  in  Figures  33a  thru  33d.  WES2DT  was  also  used  to  make  the  same  analysis. 

Figures  34a  thru  34d  show  resultant  displacement  vectors  and  Figures  35a  thru 

35d  show  principal  stress  vectors  from  the  WES2DT  analysis. 

Two-dimensional  gravity  and  thermal 
loading  including  creep  and  shrinkage 

69.  ABAQUS  was  used  to  make  this  analysis  using  FE  grid  1  and  modulus 

subroutine  UMAT2 .  Displaced-structure  plots  for  this  analysis  are  shown  in 

Figures  36a  thru  36d.  Maximum  principal  stress  contours  are  shown  in  Figures 

37a  thru  37b. 

Three -dimensional 
gravity  loading 

70.  ABAQUS  was  used  to  perform  a  3-D  gravity  turn-on  analysis  of  mono¬ 
lith  L-13  using  the  FE  grid  in  Figure  8a.  The  displaced  structure  at  the 
element  faces  in  the  5  element  planes  in  Figure  8d  are  shown  in  Figures  38a 


thru  38e. 
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Discussion  of  Results 

Gravity  turn-on  analyses 

71.  Comparisons  of  gravity  turn-on  analyses  of  monolith  L-13  made  with 
E  -  3.12  x  10^  psi  and  E  -  4.8  x  10^  psi  seen  in  Figures  16  and  17,  respec¬ 
tively,  showed  both  the  displacements  and  stresses  to  be  in  close  agreement. 
Maximum  principal  stresses  are  slightly  higher  for  E  —  4.8  x  10^  psi  as  expect¬ 
ed.  Also,  comparisons  of  displacements  between  these  2-D  runs  and  the  3-D 
gravity  turn-on  run  seen  in  Figure  38  closely  agree. 

Incremental  construction 
with  gravity  loading  only 
ABAOUS  versus  WES2DT 

72.  Incremental  construction  analyses  of  monolith  L-13  with  gravity 
loading  only  were  conducted  with  ABAQUS  and  WES2DT.  Comparison  of  displacement 
results  in  Figures  18  and  20  show  the  results  to  be  very  similar.  The  largest 
deflection  of  the  base  slab  occurs  midway  between  the  centerline  and  the  outer 
edge  of  the  monolith  after  placement  of  two  lifts  (Figures  18a  and  20a).  This 
is  because  the  loading  is  uniform  and  the  pile  support  is  weaker  in  this  area. 

As  concrete  in  the  wall  was  placed  and  the  loading  was  no  longer  uniform, 
displacements  of  the  slab  under  the  wall  increased  and  maximum  displacements  in 
the  slab  shifted  outward  toward  the  edge  of  the  slab  (Figures  18b-c  and  20b-c) . 

73.  Stress  comparisons  between  the  ABAQUS  and  WES2DT  programs  required 
comparing  principal  stress  in  contour  and  vector  plots  (Figures  19  and  21), 
respectively,  because  these  were  the  available  modes  of  displaying  principal 
stress  from  the  two  programs.  Although  comparing  data  from  contour  and  vector 
plots  was  not  an  easy  task,  it  did  appear  that  maximum  tensile  stresses 
occurred  in  the  same  locations  with  both  programs.  In  addition,  the  location 
of  these  maximum  tensile  stresses  coincided  with  locations  of  maximum  slab 
curvature,  namely  on  the  top  of  the  slab  at  the  centerline  and  at  the  bottom  of 


the  slab  directly  under  the  vertical  cantilever  section  of  the  wall.  The  two 
programs  use  different  methods  for  starting  the  stress  calculations  in  a  new 
lift.  The  slight  differences  in  stress  distributions  calculated  by  the  two 
programs  in  some  areas  can  be  attributed  to  the  different  lift  start-up  methods 
and  are  amplified  by  rapidly  increasing  values  of  E  at  early  ages. 


Gravity  turn-on  versus 
incremental  construction 

74.  Comparisons  between  whole -structure  gravity  turn-on  and  incremental 
construction  simulation  with  gravity  loading  only  with  ABAQUS  showed  predict¬ 
able  results.  In  gravity  turn-on  analyses,  the  maximum  displacements  in  the 
slab  were  nearly  equal  under  the  wall  and  lower  at  the  centerline  of  the 
monolith.  In  incremental  construction  analyses,  maximum  displacements  in  the 
slab  occurred  directly  under  the  vertical  cantilever  section  of  the  wall  with 
lower  displacements  near  the  outer  edge  of  the  slab  and  at  the  centerline.  The 
displacements  in  the  wall  were  characteristically  different  between  gravity 
turn-on  and  incremental  construction  analyses.  During  incremental  construction, 
each  successive  lift  is  placed  up  to  its  planned  elevation.  Therefore,  the 
displacement  of  the  top  of  each  lift  decreases  at  higher  elevations  in  the 
structure.  Conversely,  in  the  gravity  turn-on  analysis,  the  maximum 
displacement  occurs  at  the  top  of  the  wall  because  all  displacements  accumulate 
at  the  highest  elevation.  These  results  have  thus  exhibited  known  modes  of 


displacement  for  gravity  turn-on  and  incremental  construction  analyses. 


75.  Stress  contours  resulting  from  ABAQUS  calculations  of  gravity  turn¬ 
on  and  incremental  construction  were  very  similar  as  expected.  Stress  values 
within  the  base  slab  were  slightly  higher  due  to  use  of  a  constant,  mature  E 
(28-day  value)  whereas  the  incremental  construction  analysis  used  aging  E 
values  which  are  lower  at  early  age. 
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76.  When  thermal  effects  due  to  heat  of  hydration  of  cement  were  added 
to  the  loading,  the  predominant  changes  in  the  response  of  the  structure  were 
an  elongation  of  the  base  slab  and  a  downward  curvature  of  the  outer  end  of  the 
base  slab.  These  effects  were  observed  from  the  both  the  ABAQUS  and  WES2DT 
results.  This  response  is  illustrated  by  comparing  Figures  22,24,26,28,  and 
30.  Elongation  of  the  slab  is  due  to  thermal  expansion.  The  curvature  of  the 
base  slab  occurs  as  a  result  of  differential  thermal  expansion  between  a  new 
lift  and  previous  lifts.  The  concrete  temperature  rise  data  shown  in  Figure  4a 
shows  that  a  new  lift  can  undergo  thermal  expansion  equivalent  to  a  potential 
temperature  rise  of  35°  F  during  the  first  five  days  after  placement.  The 
previous  lift  during  the  same  time  period  (its  second  5  days)  can  only  experi¬ 
ence  thermal  expansion  equivalent  to  a  5°  F  temperature  rise.  The  differential 
longitudinal  expansion  due  to  the  temperature  rise  in  each  lift  causes  the 
downward  curvature  of  the  outer  portion  of  the  slab.  The  curvature  is  not 
actually  as  pronounced  as  the  temperature  differential  between  the  two  lifts 
would  indicate  because  the  E,  hence  stiffness,  of  the  newer  lift  is  lower  than 
that  of  the  previous  lift.  Effectively,  only  a  part  of  the  temperature  differ¬ 
ential  thus  causes  curvature  of  the  slab. 

77.  The  addition  of  thermal  loading  caused  an  increase  in  tensile 
stresses  on  the  top  surface  of  the  base  slab  near  the  outer  edge  after  place¬ 
ment  of  lift  2.  Tensile  stresses  also  increased  along  the  vertical  wall 
surfaces  and  compressive  stresses  developed  in  the  center  of  wall  masses. 

These  stresses  in  the  wall  develop  due  to  the  differential  expansion  from  the 
thermal  gradients  between  the  cooler  outer  surfaces  and  the  waim  interiors.  In 
ABAQUS  maximum  principal  stress  contour  plots  include  z-direction  stresses 
which  result  from  the  absolute  restraint  to  thermal  expansion  in  its  plane 


strain  formulation  in  two  ways.  First, 


z-direction  compressive  stresses 


increase  x-  and  y-  direction  tensile  stresses  due  to  Poisson  effect  in  ABAQUS . 
And  z-direction  out-of-plane  stresses  may  be  included  in  2-D  maximum  principal 
stress  plots.  Therefore,  z-direction  stresses  may  be  mixed  in  with  x  -  y  plane 
stresses  in  the  ABAQUS  principal  stress  plots.  It  is  understood  that  subse¬ 
quent  ABAQUS  versions  will  include  only  in-plane  stresses  in  these  plots. 

78.  Comparison  of  principal  stress  contours  (Figures  23  and  25)  using 
ABAQUS  grids  1  and  2  (Figures  7a  and  7b)  showed  virtually  identical  results. 
Comparisons  of  displacements  (Figures  22  and  24)  also  showed  very  close  agree¬ 
ment.  It  was  concluded,  based  upon  these  results,  that  the  coarse  mesh  of  FE 
grid  1  was  adequate  for  all  subsequent  analyses  of  monolith  L-13. 

79.  Comparisons  were  also  made  between  analyses  results  from  ABAQUS  and 
from  WES2DT  that  included  thermal  loading.  Comparisons  between  results  from 
the  two  programs  must  be  tempered  by  the  realization  that  ABAQUS  plane  strain 
calculations  effectively  fully  restrain  z-direction  expansion  and  WES2DT  does 
not.  Consequently,  z-direction  stresses  develop  in  ABAQUS  which  through 
Poisson  effect  modify  x-  and  y-direction  stresses.  In  U’ES2DT,  z-direction 
stresses  are  not  calculated  so  that  they  are  effectively  zero  with  no  resulting 
Poisson  effect.  In  addition,  stress  plots  in  ABAQUS  may  include  z-direction 
stresses.  With  these  points  in  mind,  examination  of  displacement  plots  from 
the  two  programs  (Figures  22  and  30)  showed  very  similar  deflection  patterns  at 
all  locations.  Comparison  of  maximum  principal  stress  plots  (Figures  23  and 
31)  showed  tensile  stress  results  that  compared  very  well  in  location  and 
magnitude.  When  maximum  principal  stresses  were  compressive,  stresses  were 
generally  lower  from  ABAQUS.  This  mav  be  the  result  of  the  different 
considerations  for  plain  strain  calculations  by  ABAQUS  and  WES2DT . 

80.  All  previously  described  analyses  with  ABAQUS  employed  the  modulus 


subroutine  UMAT1  in  which  the  E  versus  age  function  was  entered  as  tabular 
input.  ABAQUS  modulus  subroutine  UMAT2  (aging  creep  model)  enters  E  versus 
age  data  by  a  mathematical  function  representation.  When  UMAT2  was  first  used, 
several  initial  comparative  FE  analyses  were  run  to  identify  any  differences 
between  UMAT1  and  UMAT2  without  employing  creep  or  shrinkage  to  establish  bases 
of  comparisons  between  the  two  subroutines.  UMAT2  incorporates  E  as  a 
function  of  temperature  as  well  as  age.  Comparative  runs  were  made  to  identify 
the  differences  introduced  for  the  temperature  dependence  of  E  for  the  range 
of  temperatures  encountered  in  these  analyses.  The  displaced  grid  (Figures  26 
and  28)  and  principal  stress  contour  (Figures  27  and  29)  plots  show  that 
temperature  dependence  of  E  does  not  appreciably  affect  the  results.  Compar¬ 
isons  of  tha  E  versus  age  representations  used  in  UMAT1  and  UMAT2  were  made 
by  comparing  displaced  grid  (Figures  22  and  28)  and  principal  stress  contour 
(Figures  23  and  29)  plots.  These  results  show  that  elongation  and  downward 
curvature  of  the  base  slab  (lifts  1  and  2)  are  slightly  greater  for  the  UMAT2 . 
In  locations  of  maximum  tensile  stresses,  UMAT2  produced  stresses  up  to  100 
percent  higher  than  those  by  UMAT1  after  five  days  after  placement  of  lift  2. 
These  higher  stresses  from  UMAT2  were  determined  to  be  the  direct  result  of 
much  higher  values  of  E  during  the  first  day  in  the  functional  representation 
used  in  UMAT2 .  The  method  for  expressing  E  versus  time  in  the  aging  creep 
model  will  require  modification  in  the  future  to  better  express  E  values  for 
the  first  day  or  so  of  age.  This  could  not  be  done  in  time  to  be  used  in  this 
study.  The  remainder  of  the  analyses  were  made  using  ABAQUS  modulus  subroutine 
UMAT2  with  the  realization  that  early-age  stresses  would  be  higher  than  normal. 

Two - d i r ?ns i ona 1  gravity  and 
thermal  loading  including  creep 


81.  Analyses  of  L-13  which  included  gravitv  and  thermal  loading 
creep  were  conducted  using  both  ABAOUS-UMAT2  with  FF.  grid  1  and  WES.T'T 
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effects  of  creep  changed  the  deflected  response  of  the  monolith.  Although  the 
displacement  trends  (Figures  32  and  34)  produced  by  ABAQUS  and  WES2DT  compared 
remarkably  well,  magnitudes  of  deflection  were  higher  with  ABAQUS. 

82.  When  comparing  the  principal  stress  plots  without  creep  for  ABAQUS 

and  WES2DT  (Figures  23  and  31)  to  principal  stress  plots  for  the  two  programs 

with  creep  (Figures  33  and  35),  the  following  observations  were  made.  The 

relaxation  of  stresses  due  to  creep  with  ABAQUS -UMAT2  reduced  peak  principal 

stresses  in  the  base  slab  from  300  psi  (Figure  29a)  to  150  psi  (Figure  33a)  at 

five  days  after  placement  of  lift  2.  At  5  days  after  placement  of  lift  4,  peak 

principal  stresses  in  the  base  slab  were  reduced  from  250  psi  (Figure  23b)  to 

50  psi  (Figure  33b).  With  WES2DT  the  similar  comparison  showed  reduction  of 

peak  principal  stresses  due  to  creep  in  the  base  slab  from  150  psi  (Figure  31a) 

to  115  psi  (Figure  35a)  at  5  days  after  placement  of  lift  2  and  from  90  psi 

(Figure  31b)  to  50  psi  (Figure  35b)  at  5  days  after  placement  of  lift  4.  These 

comparisons  show  that  initial  stresses  are  much  higher  and  relaxation  due  to 

creep  is  much  greater  with  ABAQUS-UMAT2  than  with  WES2DT .  The  higher  initial 

stresses  are  consistent  with  the  fact  that  the  early-age  E  values  are  greater 

with  UMAT2 .  As  discussed  in  paragraph  28,  it  was  known  that  the  creep  data 

incorporated  into  UMAT2  would  produce  slightly  greater  creep  relief  than  the 

data  used  by  WES2DT.  Also,  since  the  highest  rate  of  creep  occurs  when  these 

stresses  are  high,  a  proportionately  larger  amount  of  stress  relief  resulted 

with  ABAQUS.  At  later  times  as  the  initial,  high  thermal  gradients  and  the 

resulting  thermal  stresses  are  reduced,  the  excessively  relaxed  stresses  are 

further  reduced.  This  process  accounts  for  the  large  decrease  in  stresses  seen 

in  the  base  slab  in  Figure  35b. 

Two-dimensional  gravity  and  thermal 
loading  including  creep  and  shrinkage 


83.  An  ABAQUS  analysis  was  conducted  on  L-13  using  UMAT2  with  gravity 


and  thermal  loading  including  creep  and  shrinkage.  The  addition  of  this 
autogenous  shrinkage  caused  noticeably  different  results.  Comparison  of  ABAQUS 
analyses  with  creep,  but  with  and  without  shrinkage  showed  both  different 
displacement  and  principal  stress  response.  Comparison  of  displaced  grids 
(Figures  32  and  36)  for  analyses  with  and  without  shrinkage,  respectively, 
shows  that  the  addition  of  shrinkage  cancels  some  of  the  thermal  elongation  of 
lift  2  and  reduces  the  differential  thermal  expansion  between  lifts  1  and  2 
that  had  produced  downward  curvature  of  the  outer  end  of  the  base  slab. 
Subsequently,  further  shrinkage  caused  deflections  of  the  slab  directly  under 
the  verticall  cantilever  section  of  the  wall  to  increase  substantially  (Figures 
36a  -  36d) .  Comparison  of  maximum  principal  stresses  showed  that  at  5  days 
after  placement  of  lift  2,  stresses  in  the  base  slab  were  greater  and  nearly 
symmetrical  about  the  lift  1-lift  2  interface  when  shrinkage  was  applied.  In 
fact,  stresses  increased  in  several  areas  of  the  wall.  Although  the  creep  and 
shrinkage  data  used  in  these  ABAQUS  FE  analyses  are  not  based  upon  the  results 
of  tests  of  Lock  and  Dam  26R  concrete  which  may  cause  the  results  to  deviate 
slightly,  it  is  evident  from  these  analyses  that  inclusion  of  creep  and 
shrinkage  are  necessary  for  proper  thermal  stress  evaluation. 

Pile  loads 

84.  The  load  in  a  pile  is  directly  proportional  to  the  deflection  of 
the  pile  head,  therefore,  the  distribution  of  pile  loads  underneath  the  base 
slab  is  dependent  principally  upon  the  base  slab  deflections.  Table  6  gives 
the  vertical  pile  loads  from  analyses  using  ABAQUS  and  WES2DT.  Also  included 
are  pile  loads  supplied  by  St.  Louis  District.  The  ABAQUS  results  using  a 
Young's  modulus  modified  by  ACI  agree  very  closely  with  the  St.  Louis  District 
values  in  a  gravity  turn-on  analysis.  Since  the  ABAQUS  and  WES2DT  modulus 
values  were  not  to  be  modified  for  the  remaining  analyses,  an  ABAQUS  analysis 
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was  also  made  with  the  unmodified  28-day  Young's  modulus.  The  stiffer  concrete 
slightly  redistributed  the  pile  loads  with  the  higher  loads  still  being  located 
under  the  wall  at  the  outer  edge  of  the  base  slab.  The  maximum  pile  loads  from 
these  analyses  were  all  less  than  200  kips/pile. 

85.  When  the  wall  was  constructed  incrementally  with  only  gravity  loads 
applied,  the  ABAQUS  and  WES2DT  programs  predicted  more  significant 
redistribution  of  the  pile  loads  with  maximum  values  of  200-210  kips/pile 
occurring  underneath  the  wall. 

86.  As  thermal  loading  was  added,  the  ABAQUS -UMAT1  results  for  grid  1 
and  grid  2  were  almost  identical.  The  results  of  ABAQUS  using  the  new  modulus 
subroutine  UMAT2 ,  for  both  modulus  as  a  function  of  temperature  and  constant 
with  temperature,  gave  very  similar  results.  No  further  analyses  were  made 
with  the  modulus  held  constant  with  temperature  since  it  caused  no  significant 
change  in  pile  loads.  However  the  predicted  pile  loads  using  UMAT2  were  redis¬ 
tributed  from  the  ones  using  UMAT1 .  The  larger  pile  loads  were  reduced  under 
the  wall  and  moved  farther  toward  the  outer  edge  of  the  base  slab.  Maximum 
pile  loads  were  approximately  230  kips/pile  for  the  ABAQUS -UMAT2  analysis. 
WES2DT  resuLts  showed  close  agreement  in  magnitude  and  trend  with  the  ABAQUS 
results  using  UMAT1 

87.  When  creep  was  included  in  the  analysis,  both  ABAQUS  and  WES2DT 
redistributed  the  pile  loads.  Both  gave  lower  pile  loads  at  the  centerline, 
lower  loads  at  the  outer  edge,  and  higher  loads  under  the  wall.  The  ABAQUS 
results  showed  more  extreme  changes  in  each  instant  WES2DT  results  ranged 
from  66  kips/pile  at  the  centerline  to  213  kips/pile  under  the  wall.  ABAQUS 
results  showed  Ub  kips/pile  at  the  centerline  and  258  kips/pile  under  the  wall. 
These  ABAQUS  values  are  10  -  15  percent  higher  than  WES2DT  values.  The  higher 
pile  loads  can  be  attributed  to  higher  creep  occurring  in  ABAQUS  at  earlv  times 
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after  placement  of  a  lift  when  the  modulus  is  excessively  high.  The  VES2DT 
pile  loads  are  very  similar  in  magnitude  and  distribution  to  those  obtained 
earlier  when  only  gravity  loading  was  included  in  an  incremental  construction 
analysis . 

88.  When  shrinkage  was  added  to  the  ABAQUS  analysis,  the  pile  loads 
were  slightly  redistributed  from  those  where  only  creep  was  included.  Higher 
values  were  predicted  under  the  vertical  cantilever  section  of  the  wall  and 
lower  values  at  the  outer  edge  of  the  base  slab. 

89.  While  the  gravity  turn-on  analysis  method  gives  adequate  pile  load 
input  for  use  in  design,  the  incremental  construction  method  should  give  loads 
closer  to  the  actual  values  since  it  more  closely  models  the  actual  construc¬ 
tion  sequence.  However,  the  incremental  method  will  be  slightly  more  expensive 
and  difficult  to  perform  since  data  to  define  the  modulus  as  a  function  of  time 
is  needed  and  a  series  of  solutions  is  performed. 

Cracking  potential 

90.  Comparisons  of  results  were  made  to  evaluate  the  potential  for 
cracking  of  the  concrete  in  the  analyses  that  were  made  in  this  investigation. 
Modulus  of  rupture  test  results  from  Lock  and  Dam  26R,  Phase  I  mixture  4c 
concrete  yielded  values  of  124  psi,  280  psi,  300  psi,  and  464  psi  at  1-,  3-,  7- 
,  and  28-days  agt»,  respectively.  If  these  values  are  used  as  a  simple  cracking 
initiation  limit  for  tensile  stresses,  the  following  conclusions  can  be  made 
regarding  crack  potential  in  monolith  L-13. 

91.  First,  maximum  principal  stresses  computed  in  WES2DT  analyses  did 
not  exceed  150  psi  in  tension  under  any  state  of  loading  ar.d  this  level  was 
reached  after  3-days  age  without  any  stress  relief  due  to  creep.  With  creep 
peak  tensile  stresses  were  20-40  percent  less.  Consequently,  maximum  tensile 


stresses  were  less  than  50  percent  of  modulus  of  rupture  at  anv  age  even 


without  benefit  of  creep. 

92.  Secondly,  maximum  principal  tensile  stresses  computed  in  ABAQUS 
analyses  using  UMAT1,  which  employed  the  identical  modulus  versus  time 
relationship  as  WES2DT  analyses,  did  not  exceed  200  psi  and  only  after  3-days 
age.  Even  without  the  benefit  of  stress  relief  through  creep,  which  should 
relax  stresses  by  20-40  percent,  these  peak  tensile  stresses  are  only  around  70 
percent  of  modulus  of  rupture  at  3 -day. 

93.  Finally,  maximum  principal  tensile  stresses  computed  in  ABAQUS 
analyses  using  UMAT2 ,  which  contained  an  excessively  high  modulus  function  for 
the  first  day  or  so  of  age,  exceeded  300  psi  at  3-  to  5-days  age  without  creep 
and  150  psi  with  creep.  Considering  that  the  computed  early-age  stresses  were 
high  due  to  the  abnormally  high  initial  modulus  used  in  UMAT2 ,  it  is  probable 
that  values  of  100  psi  or  less  with  creep  at  3-days  age  can  be  expected.  This 
reinforces  the  conclusion  reached  with  the  WES2DT  analyses  that  peak  tensile 
stresses  are  less  than  50  percent  of  modulus  of  rupture  in  this  structure. 
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Figure  16a. 


Displaced  structure  from  gravity  turn-on  analysis 

E  -  3.12  x  106  psi 
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Figure  16b.  Maximum  principal  stress  contours  for  gravity  turn-on  analy 
E  -  3.12  x  106  psi 
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Figure  17a.  Displaced  structure  from  gravity  turn-on  analysis, 
E  -  4 . 8  x  10  psi 
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Figure  18b.  Displaced  structure  5  davs  after  lift  4  is  placed 
gravity  loading  o:ilv,  no  creep,  using  program  ABAQUS 


with  first  grid 
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Figure  18c.  Displaced  structure  5  days  after  lift  5  is  placed 
gravity  loading  only,  no  creep,  using  program  AilAQUS 

with  first  grid 
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Figure  18d.  Maximum  principal  stress  (tensile)  contour 
days  after  placement  of  lift  9,  gravity  loading  on 
using  program  ABAQllS  (UMAT  1)  first  grid 


Figure  19a.  Maximum  principal  stress  (tensile)  contours  in  structure  5 
days  after  placement  of  lift  2,  gravity  loading  only,  no  creep,  using 

program  ABAQUS  with  first  grid 
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Figure  19d.  Displaced  structure  7  days  after  lift  9  is  placed,  gravity 
loading  only,  no  creep,  using  program  ABAQUS  (UMAX  1)  with  first  grid 
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Figure  22a.  Displaced  structure  5  days  after  lift  2  is  placed,  thermal  and  grav 
loading,  no  creep,  using  program  ABAQUS  witn  first  grid  (L»1A 
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Figure  22b.  Displaced  structure  5  days  after  lift  U  is  placed,  thermal  and  gra 
loading,  no  creep,  using  program  ABAQUS  wit:i  first  grid  (UMAT  1) 
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Figure  22c. 
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loading,  no  creep,  using  program  ABAQUS  with  firs 


thermal  and  gravity 
grid  (UMAT  1) 
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Figure  23a.  Maximum  principal  stress  (tensile)  contours  in  structure  3  days  after 
placement  of  lift  2,  thermal  and  gravity  loading,  no  creep,  using  program  ABAQL'S 

with  first  grid 
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Figure  25c.  Maximum  principal  stress  (tensile)  contours  in  structure  5  days  after 
placement  of  lift  5,  thermal  and  gravity  loading,  no  creep,  using  program  ABAQl’S 

with  second  grid 
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Figure  26a.  Displaced  structure  5  days  after  placement  of  lift  2,  E  constant 
with  temperature,  using  ABAQUS  (UMAT2) . 
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Figure  26c.  Displaced  structure  5  days  after  placement  of  lift  5, 
with  temperature,  using  ABAQUS  (UMAT2) . 
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Figure  27a.  Maximum  principal  stress  contours  5  days  after  placement  o 
lift  2,  E  constant  with  temperature,  using  ABAQUS  (UMAT2). 
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Figure  27b.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  4,  E  constant  with  temperature,  using  ABAQUS  (UMAT2) . 
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Figure  27c.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  5,  E  constant  with  temperature,  using  ABAQUS  (UMAT2) . 
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Figure  27d.  Maximum  principal  stress  contours  7  days  after  placement  of 
lift  9,  E  constant  with  temperature,  using  ABAQUS  (UMAT2). 
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Displaced  structure  5  days  after  placement  of  lift  2,  E  is 
function  of  temperature,  using  ABAQUS  (UMAT2) . 
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Figure  28b.  Displaced  structure  5  days  after  placement  of  lift  4,  E  is  a 
function  of  temperature,  using  ABAQUS  (UMAT2). 
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Figure  28c.  Displaced  structure  5  days  after  placement  of  lift  5,  E  is  a 
function  of  temperature,  using  ABAQUS  (UMAT2). 


i.I"?  -  m  1  3  -  COARSE  MODEL  -NE'J  UMAT.NO  SHR  I NK,  NOCREEP 

STEP  3S  [HCPEPEHT  2  flanOUS  yE»5tO»  4-S-I17 


Figure  28d.  Displaced  structure  7  days  after  placement  of  lift  9,  E  is 
function  of  temperature,  using  ABAQl'S  (UMAT2). 
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Figure  29a.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  2,  E  is  a  function  of  temperature,  using  ABAQUS  (UMAT2) 
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Figure  29b.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  4,  E  is  a  function  of  temperature,  using  ABAQl'S  (UMAT2) 
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Figure  29d.  Maximum  principal  stress  contours  7  days  after  placement  of 
lift  9,  E  is  a  function  of  temperature,  using  ABAQUS  (UMAT2) 
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Figure  32b.  Displaced  structure  5  days  after  placement  of  lift  U  including 
creep,  using  ABAQUS . 
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lift  5  is  placed,  thermal  and  gravity  loading 
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Figure  37a.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  2  including  creep  and  shrinkage,  using  ABAQUS  (UMAT2). 
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Maximum  principal  stress  contours  5  days  after  placement  of 
lift  A  including  creep  and  shrinkage,  using  ABAQUS  (UMAT2). 
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Figure  37c.  Maximum  principal  stress  contours  3  days  after  placement  of 
lift  5  including  creep  and  shrinkage,  using  ABAQUS  ( UMAT 2 ) . 
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Figure  3Sa.  Displaced  structure  for  3-dimensional  analysis  of  Monolith  13 
at  node  plane  1  (centerline). 
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Figure  33b.  Displaced  structure  for  3 - dimens iorial  analysis  of  Monolith  13 
at  node  plane  2. 
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PART  IX:  THERMAL  ANALYSES  OF  MONOLITH  L-17 


Presentation  of  Results 


Two-dimensional 
temperature  analysis 

94.  Two-dimensional  temperature  analyses  of  monolith  L-17  were  also 

performed  using  both  the  ABAQUS  and  WES2DT  programs.  The  2-D  FE  models  used  by 

ABAQUS  (Figure  7d)  and  WES2DT  (Figure  7e)  were  based  upon  a  cross-section  taken 

at  Station  26  +  13  which  contained  the  smallest  amount  of  internal  voids. 

This  cross-section  was  expected  to  experience  the  largest  temperature  rise  in 

the  monolith.  The  WES2DT  grid  only  modeled  the  structure  through  lift  13. 

During  the  temperature  analysis,  the  top  surface  of  lift  13  was  insulated  after 

the  time  of  placement  of  lift  14  to  simulate  the  effects  of  the  additional 

concrete.  The  gravity  loading  effects  past  lift  14  were  modeled  with  pressure 

loads  applied  to  the  top  of  lift  13  to  simulate  placement  of  lifts  14  through 

16  during  the  incremental  construction  stress  analysis.  Figures  39a  thru  39e 

show  results  from  the  ABAQUS  program  at  five  stages  of  construction.  Results 

from  the  WES2DT  program  are  shown  in  Figures  40a  thru  40d  for  the  same  first 

four  of  five  stages  of  construction  used  for  presenting  ABAQUS  data.  The 

presentation  of  WES2DT  temperature  results  beyond  the  time  of  placing  lift  13 

was  omitted. 

Three  dimensional 
temperature  analysis 

95.  Although  a  3-D  FE  model  of  L-17  (Figure  7e)  for  use  with  the  ABAQUS 
program  was  completed,  the  analysis  could  not  be  run.  The  temperature  analysis 
was  not  made  due  to  problems  with  the  ABAQUS  program  in  handling  the  required 
number  of  element  sets  and  the  excessive  costs  of  making  the  computer  runs. 


190 


Discussion  of  Results 


Two-dimensional 
temperature  analysis 

96.  The  results  from  ABAQUS  and  WES2DT  show  very  good  agreement.  The 
temperature  contours,  Figures  39  and  40,  give  the  same  trends  and  are  in  the 
identical  locations  for  times  through  placement  of  lift  13.  Slightly  larger 
areas  are  shown  for  the  100°  F  contours  from  ABAQUS  in  Figures  39b  and  39d  that 
from  UES2DT  in  Figures  40b  and  40d.  The  ABAQUS  analysis  shows  the  maximum 
temperature  of  105°  -  110°  F  occurs  after  placement  of  lift  16  in  the  center  of 
the  massive  wall  seccion.  Since  the  WES2DT  model  only  extended  through  lift  13 
the  corresponding  value  is  not  given. 


191 


1.0.  VALUE 

1  «8. S0E»01 

2  •T.30E*01 

3  .7.S0E»01 
it  .a.30E*oi 
s  »a.50E»01 

8  ♦9.39£»01 
7  *9.50E*01 

9  *1 . 30£*02 
9  • ! . 35E*02 

10  ♦  !. 1 0£*02 
11  •!  15E*02 


L IFT4  -  M17 

STEP  B  INCREHENT  S 


20  MODEL 

ABAQUS  VERSION  N-S-U7 


Figure  39a. 


Temperature  contours  calculated  using  ABAQUS  in  Monolith  17  at 
5  days  after  placement  of  lift  4. 
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PART  X:  STRESS  ANALYSES  OF  MONOLITH  L-17 


97.  The  ABAQUS  program  was  used  to  perform  the  2-D  stress  analysis  of 


monolith  L-17  for  both  instantaneous  gravity  turn-on  of  the  entire  structure 


and  an  incremental  construction  build-up  sequence  of  the  structure.  Figures 


41a  and  41b  show  the  displaced  structure  and  maximum  principal  (tensile)  stress 


contours,  respectively,  for  the  gravity  turn-on  analysis  for  a  Young's  Modulus 


(E)  of  3.12  x  10°  psi.  Figures  42a  and  42b  give  the  same  results  for  E  -  4.80 


x  10°  psi.  These  modulus  values  correspond  to  the  ACI  modified  and  normal 


compressive  strength  of  the  concrete  at  28  days,  respectively.  Figures  43a 


thru  43c  and  44a  thru  44c  show  the  displaced  structure  and  maximum  principal 


stress  contours,  respectively,  for  three  stages  of  construction  using  incremen¬ 


tal  construction  sequencing.  This  incremental  construction  analysis  uses  the 


material  modulus  calculation  routine  UMAT1 . 


98.  WES2DT  was  used  to  make  an  incremental  construction  analysis.  As 


described  in  paragraph  94,  a  FT  model  was  prepared  through  lift  13.  Gravity 


loads  for  lifts  14  through  17  were  simulated  by  applying  appropriate  surface 


pressures  to  the  top  of  lift  13.  Figures  45a  thru  45c  and  46a  thru  46c  show 


the  displacement  vector  plots  and  principal  stress  vector  plots,  respectively, 


from  this  analysis.  Figures  45c  and  46c  represent  the  conditions  at  5  davs 


after  lift  16  is  placed  which  corresponds  to  ABAQUS  analyses  output  in  Figures 


43c  and  44c.  All  subsequent  plots  of  WES2DT  results  will  be  similarly  shown. 


Two  -  dimens iona 1 

gravity  and  thermal  loan  in; 


99.  The  ABAQUS  program  was  used  to  ..erf'".':;;  t  he  2  -D  stress  analysis 


including  gravity  and  temperature  loading  us  i  tv,  -he  modulus  routine  U.MA71  .  Th 


■  Ail  A  JiV  m  Aa  Vm  ■_«  t  a  ",  Aa  .Aa  Aj  A_m  ~ 


displaced  structure  plots  are  given  in  Figures  47a  thru  47c.  The  maximum 
principal  stress  contours  are  given  in  Figures  48a  thru  48c.  The  modulus 
subroutine  UMAT2  (concrete  aging  creep  model)  was  used  with  ABAQUS  to  make  an 
analysis  with  gravity  and  thermal  loading.  The  displaced  structure  and  maximum 
principal  stress  contours  are  given  in  Figures  49a  thru  49c  and  50a  thru  50c, 
respectively.  A  WES2DT  analysis  was  not  made  for  this  same  gravity  and  tempe¬ 
rature  loading  only. 

Two-dimensional  gravity  and 
thermal  loading  including  creep 

100.  ABAQUS  was  used  in  incremental  construction  analyses  of  the  struc¬ 
ture  with  the  modulus  subroutine  UMAT2  including  creep.  Displaced- structure 
plots  are  shown  in  Figures  51a  thru  51c.  Maximum  principal  stress  contours  are 
shown  in  Figures  52a  thru  52c.  WES2DT  was  also  used  to  make  the  same  analysis. 
Figures  53a  thru  53c  show  resultant  displacement  vectors  and  Figures  54a  thru 
54c  show  principal  stress  vectors  from  the  WES2DT  analysis. 

Two-dimensional  gravity  and  thermal 


loading  including  creep  and  shrinkage 


101.  ABAQUS  was  used  to  make  this  analysis  also  with  modulus  subroutine 
UMAT2 .  Displaced- structure  plots  for  this  analysis  are  shown  in  Figures  55a 
thru  55c.  Maximum  principal  stress  contours  are  shown  in  Figures  56a  thru  56c. 


Discussion  of  Results 


Gravity  turn-on  analyses 

102.  Comparisons  of  gravity  turn-on  analyses  of  monolith  L-17  made  with 
E  ”  3.12  x  106  psi  and  E  -  4.8  x  10^  psi  seen  in  Figures  41and  42,  respec¬ 
tively,  showed  both  the  displacements  and  stresses  to  be  in  close  agreement. 
Maximum  principal  stresses  are  slightly  higher  for  E  —  4.8  x  10^  psi  as  expect¬ 


ed  especially  in  the  top,  center  of  the  base  slab. 


103.  Incremental  construction  analyses  of  monolith  L-17  with  gravity 
loading  only  were  conducted  with  ABAQUS  and  WES2DT.  Comparison  of  displacement 
results  in  Figures  43  and  45  show  the  results  to  be  very  similar.  The  largest 
deflections  in  the  base  slab  occur  all  along  its  entire  length  except  near  the 
center  of  the  wall  zone  after  placement  of  the  first  four  lifts  (Figures  43a 
and  45a) .  This  is  because  the  loading  is  uniform  and  the  pile  support  is 
stronger  under  the  wall.  As  concrete  in  the  lower  wall  region  was  placed, 
deflections  in  the  base  slab  became  more  uniform  (Figures  43b  and  45b) .  As 
wall  placement  continued,  displacements  of  the  slab  under  the  wall  increased 
and  maximum  displacements  in  the  slab  shifted  outward  toward  the  edge  of  the 
monolith  (Figures  43c  and  45c). 

104.  Stress  comparisons  between  the  ABAQUS  and  WES2DT  programs  required 
comparing  principal  stress  in  contour  and  vector  plots  (Figures  44  and  46) , 
respectively,  as  was  described  earlier  in  discussion  of  L-13.  Maximum  tensile 
stresses  were  generally  in  the  same  locations  of  the  models  with  both  programs. 


105.  Comparisons  between  whole -structure  gravity  turn-on  and  incremen¬ 
tal  construction  simulation  with  gravity  loading  only  with  ABAQUS  showed  pre¬ 
dictable  results.  In  both  analyses,  the  maximum  displacements  of  the  base 
occurred  after  the  simulation  reached  full  height.  Displacements  were  less  at 
the  centerline  of  the  slab  and  greatest  near  the  outer  edge  ofthe  monolith 
under  the  wall.  Actually,  displacements  at  the  outer  edge  were  greater  in  the 
incremental  construction  analysis.  As  was  demonstrated  in  the  L-13  analyses, 
the  displacements  within  the  base  slab  and  wall  were  characteristically  diffe¬ 
rent  in  gravity  turn-on  and  incremental  construction  analyses.  During  increme- 
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ntaL  construction,  each  successive  lift  is  placed  up  to  its  planned  elevation. 
Therefore,  the  displacement  of  the  top  of  each  lift  decreases  at  higher  eleva¬ 
tions  in  the  structure.  Conversely,  in  the  gravity  turn-on  analysis,  the 
maximum  displacement  occurs  at  the  top  of  the  wall  because  all  displacements 
accumulate  at  the  highest  elevation.  These  results  have  thus  exhibited  known 
modes  of  displacement  for  gravity  turn-on  and  incremental  construction  ana¬ 
lyses  . 

106.  Stress  contours  resulting  from  ABAQUS  calculations  of  gravity 
turn-on  and  incremental  construction  were  very  similar  as  expected.  Stress 
values  within  the  base  slab  were  slightly  higher  for  the  gravity  turn-on  analy¬ 
sis  which  used  E  -  4.8  x  10^  psi  and  for  incremental  construction  analysis 
because  E  values  were  more  similar  while  the  wall  was  placed. 

Two-dimensional 

eravitv  and  thermal  loading 

107.  When  thermal  effects  due  to  heat  of  hydration  of  cement  were  added 
to  the  loading,  the  predominant  changes  in  the  response  of  the  structure  were, 
as  with  L-13  analyses,  an  elongation  of  the  base  slab  and  a  d<  toward  curvature 
of  >_he  outer  end  of  the  base  slab.  However,  in  L-17  the  curvature  phenomenon 
continued  into  the  wall.  These  effects  were  observed  from  the  both  the  ABAQUS 
and  WES2DT  results.  This  response  is  illustrated  by  comparing  Figures  47,  49, 
51,  and  53.  Elongation  of  the  slab  is  due  to  thermal  expansion  relative  to  the 
centerline.  As  described  in  paragraph  76,  the  curvature  occurs  as  a  result  of 
differential  thermal  expansion  between  a  new  lift  and  previous  lifts.  In  FE 
analyses  of  monolith  L-13,  the  differential  expansion  phenomenon  only  occurred 
in  the  two  lifts  of  the  base  slab.  None  of  the  higher  lifts  was  sufficiently 
long  for  any  observable  curvature  to  result.  In  monolith  L-17,  however,  there 
are  four  full  and  one  partial  full-width  lifts  in  the  base  and  nine  full-width 


lifts  in  the  wall.  The  differential  thermal  expansion  of  all  these  lifts 


contributes  to  a  curvature  in  both  the  base  slab  and  wall. 

108 .  The  addition  of  thermal  loading  caused  an  increase  in  maximum 
principal  (tensile)  stresses  along  the  bottom  of  the  base  slab  near  the  center- 
line  and  on  the  top  surface  of  the  base  slab  near  the  outer  edge  through 
placement  of  lift  4.  The  maximum  principal  stress  response  was  similar  up  to 
this  stage  of  simulated  construction  from  use  of  both  modulus  subroutines  UMAT1 
and  UMAT2  (Figures  48a  and  50a) .  The  increases  in  tensile  stresses  due  to 
thermal  loading  in  the  areas  cited  above  were  more  than  300  psi  and  100  psi, 
respectively.  Results  from  use  of  UMAT1  and  UMAT2  began  to  differ  by  the  time 
lift  16  in  the  wall  was  placed.  With  UMAT1,  the  zone  of  maximum  tensile  stress 
continued  to  be  located  along  the  bottom  of  the  base  slab  near  the  centerline. 
With  UMAT2 ,  however,  the  zone  of  maximum  tensile  stress  shifted  to  the  top  of 
the  base  slab  near  the  centerline.  This  difference  in  maximum  principal  stress 
response  is  again  attributed  to  the  higher  values  of  E  within  the  first  day 
with  UMAT2 . 

109.  The  addition  of  thermal  loading  also  increased  the  level  of  maximum 
principal  (tensile)  stresses  in  the  wall.  Tensile  stresses  increased  from  near 
zero  with  gravity  loading  only  to  about  200  psi  on  the  top  surface  of  lift  8 
with  both  UMAT1  and  UMAT2 .  As  lifts  9  through  16  were  added,  maximum  principal 
stresses  in  the  upper,  massive  portion  of  the  wall  appeared  to  coincide  with 
lift  interfaces  in  a  vertically,  alternating  pattern  of  higher  and  lower  stress 
contours  (Figures  48c  and  50c).  It  is  believed  that  the  higher  tensile  stress 
levels  that  exist  on  the  top  surface  of  each  new  lift  are  "locked  in"  as  a 
result  of  placement  of  the  successive  lift.  Lower  levels  of  maximum  tensile 
stresses  in  these  areas  were  generated  in  the  UMAT2  analysis.  This  situation 
may  be  attributed  to  the  higher,  early-age  E  values  used  by  the  subroutine 
UMAT2 .  It  is  believed  that  different  stress  fields  are  developed  when  the 
values  of  E  change  during  the  same  time  that  a  temperature  change  is  ongoing 


in  the  concrete.  The  higher  initial  E  in  UMAT2  prevents  this  same  reaction. 
A  comparable  WES2DT  analysis  with  gravity  and  thermal  loading,  but  without 
creep  was  not  conducted. 


Two-dimensional  gravity  and 
thermal  loading  including  creep 

110.  Analyses  of  L-17  which  included  gravity  and  thermal  loading  with 
creep  were  conducted  using  both  ABAQUS-UMAT2  and  WES2DT.  The  effects  of 
creep,  as  on  monolith  L-13,  changed  the  displaced  response  of  the  monolith. 

When  ABAQUS  results  using  UMAT2  with  creep  were  compared  with  those  without 
creep,  it  was  seen  that  creep  decreases  the  curvature  of  the  outer,  bottom  of 
the  base  (Figures  49  and  51).  Outward  expansion  of  the  base  slab  and  wall 
increased  with  creep,  the  direct  result  of  displaced  relaxation  of  the  lower 
lifts  to  the  elongating  force  supplied  by  the  thermally  expanding  subsequent 
lifts.  It  seems  reasonable  certain  that  this  effect  is  exaggerated  due  to  the 
higher,  early-age  E  values  used  by  UMAT2 .  With  creep,  the  concrete  appeared 
to  displace  in  an  exaggerated  manner  in  certain  locations.  This  was  the  case 
for  the  increased  vertical  displacement  of  the  outer  edge  of  the  floor  slab 
(Figures  49c  and  51c). 

111.  The  displacement  trends  produced  by  ABAQUS  and  WES2DT  did  not 
compare  (Figures  51  and  53)  as  well  as  was  seen  earlier  with  results  from  L-13. 
Again  this  may  be  partially  due  to  the  different  plane  strain  formulations 
used  in  the  two  programs  that  was  discussed  earlier.  The  larger  mass  of  this 
monolith  increased  internal  temperatures.  Therefore,  increased  potential 
existed  for  differential  expansion.  It  has  been  seen  that  this  differential 
caused  generalized  curvature  of  the  monolith.  However,  the  combined  effects  of 
high  early-age  E  ,  higher  creep  potential  than  used  by  WES2DT.  increased 
member  stiffness  due  to  size,  and  higher  temperature  rise  caused  larger  lateral 


expansion.  Actually,  the  displacements  from  the  WES2DT  analysis  more  closely 
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matched  the  ABAQUS-UMAT2  analysis  without  creep.  This  may  indicate  that  the 
amount  of  creep  in  the  ABAQUS-UMAT2  analysis  is  excessive. 

112.  When  comparing  the  ABAQUS-UMAT2  maximum  principal  stress  plots 

without  creep  (Figures  50)  to  maximum  principal  stress  plots  with  creep  (Figure 

52),  the  following  observations  were  made.  The  relaxation  of  stresses  due  to 

creep  with  ABAQUS-UMAT2  reduced  peak  principal  stresses  in  the  base  slab  from 

300  psi  (Figure  50a)  to  100  psi  (Figure  52a)  at  five  days  after  placement  of 

lift  4.  At  5  days  after  placement  of  lift  7.  peak  principal  stresses  in  the 

base  slab  were  reduced  from  around  300  psi  (Figure  50b)  without  creep  to  50  psi 

with  creep  (Figure  52b) .  Maximum  principal  stresses  in  the  vicinity  of  lift  13 

at  5  days  after  placement  of  lift  16  were  100  psi  without  creep  and  near  zero 

psi  with  creep.  In  contrast,  maximum  principal  stresses  in  the  base  slab  from 

WES2DT  with  creep  were  about  225  psi  at  5  days  after  placement  of  lift  4,  and 

approximately  175  psi  at  5  days  after  placement  of  lift  7.  Although  WES2DT  was 

not  run  without  creep  for  L-17,  the  very  large  relaxation  of  stresses  due  to 

creep  with  ABAQUS-UMAT2  and  the  higher  stresses  with  creep  from  WES2DT  act  to 

confirm  the  apparent  excessive  creep  from  UMAT2 .  This  is  the  same  conclusion 

as  reached  in  paragraph  82  for  the  analyses  of  L-13.  It  is  expected  that 

incorporating  into  UMAT2  a  corrected  E  versus  time  modulus  model  that  better 

represents  the  1-  to  2 -day  values  and  comparable  creep  data  as  were  used  in  the 

WES2DT  analysis,  lower  initial  stress  levels  and  less  creep  relaxation  would 

result  from  ABAQUS  analyses. 

Two-dimensional  gravity  and  thermal 
loading  including  creep  and  shrinkage 

113.  An  ABAQUS  analysis  was  conducted  on  L-17  using  UMAT2  with  gravity 
and  thermal  loading  including  creep  and  shrinkage.  The  addition  of  this  auto¬ 
genous  shrinkage  caused  similar,  noticeablv  different  results  as  was  seem 
earlier  in  L-13.  Comparison  of  ABAQUS  analyse 


107 


s  with  creep,  but  with  and  with- 


out  shrinkage  showed  both  different  displacement  and  principal  stress  response. 
Comparison  of  displaced  grids  (Figures  51  and  55)  for  analyses  with  and  without 
shrinkage,  respectively,  shows  that  the  addition  of  shrinkage  greatly  reduced 
the  thermal  elongation  of  successive  lifts  in  the  base  slab  and  reduced  the 
differential  thermal  expansion  between  these  lifts.  The  result  was  greatly 
reduced  downward  curvature  of  the  outer  end  of  the  base  slab  and  a  general 
shortening  of  the  slab.  Shrinkage  did  cause  the  wall  to  attain  a  curvature 
vertically  along  both  the  inner  and  outer  surfaces  (Figure  55c). 

114.  Comparison  of  maximum  principal  stresses  showed  that  at  5  days 
after  placement  of  lift  4,  peak  values  were  nearly  the  same  as  without  shrink¬ 
age.  However,  the  location  of  maximum  stresses  changed  from  the  bottom  of  the 
base  slab  without  shrinkage  to  the  top  of  the  slab  with  shrinkage  applied. 
Maximum  principal  stresses  also  increased  in  the  wall,  especially  at  times 
shortly  after  placement.  This  is  seen  in  comparisons  of  maximum  principal 
stresses  in  Figures  52b  -  52c  and  56b  -  56c  .  Stresses  in  the  last  two  lifts 
placed  with  shrinkage  are  clearly  higher.  Shrinkage  provides  an  additive 
effect  to  thermal  stresses  that  result  from  internal  restraint. 

Pile  loads 

115.  As  stated  earlier,  the  distribution  of  the  pile  loads  along  the 
bottom  of  the  base  slab  is  principally  dependent  upon  the  base  slab  deflec¬ 
tions.  Table  9  gives  the  vertical  pile  loads  from  the  different  ABAQUS  and 
WES2DT  analyses  and  those  supplied  by  the  St.  Louis  District. 

116.  The  ABAQUS  results  for  a  gravity  turn-on  analysis  using  an  ACI- 
modified  modulus  value  is  in  close  agreement  with  the  St.  Louis  District 
results.  The  discrepancies  can  he  attributed  to  minor  cross-sectional 
differences  in  the  two  cases  An  ARAQUS  anal /sis  was  also  made  using  the 


unmodified  modulus  value  since  the  remaining  incremental  analysis  would  also 


use  unmodified  values.  As  in  the  L-13  runs,  this  stiffer  concrete  slightly 
redistributed  the  pile  loads  with  the  higher  values  being  located  under  the 
wall.  Maximum  loads  were  approximately  350  kips/pile. 

117.  When  the  lock  was  constructed  incrementally  using  ABAQUS  and  WES2DT 
with  only  gravity  loading,  a  slight  redistribution  of  pile  loads  took  place. 
Smaller  loads  were  located  under  the  wall  and  larger  loads  were  located  under 
the  chamber  section  of  the  slab.  These  shifts  in  loads  were  less  than  20 
kips/pile  with  very  good  agreement  between  ABAQUS  and  WES2DT.  The  maximum  pile 
loads  occurred  at  the  same  location  with  only  a  slight  decrease  in  magnitude 
when  compared  with  gravity  turn-on  analyses. 

118.  As  thermal  loading  was  added,  the  pile  loads  shifted  more  to  the 
wall  area  with  both  UMAT1  and  UMAT2  modulus  subroutines.  With  UMAT1  the  maxi¬ 
mum  pile  load  was  460  kips/pile  while  with  UMAT2  it  was  520  kips/pile.  These 
maximum  loads  occurred  at  the  second  pile  in  from  the  outer  edge  of  the  mono¬ 
lith  in  both  cases. 

119.  When  creep  was  considered  in  the  analysis,  the  ABAQUS  and  WES2DT 
programs  gave  similar  results  except  for  piles  10  -  15.  This  is  principally 
due  to  the  WES2DT  analysis  using  a  smeared  pile  stiffness  over  a  given  are 
while  ABAQUS  used  discrete  piles.  ABAQUS  and  WES2DT  produced  larger  pile  loads 
under  the  inner  wall  area  while  having  lower  values  under  the  wall  toward  the 
outer  edge.  The  maximum  values  from  ABAQUS  was  376  kips/pile  and  again  occur¬ 
red  at  the  second  pile  from  the  outer  edge  of  the  monolith.  The  WES2DT  results 
gave  the  maximum  load  in  pile  14  under  the  inner  wall.  However,  this  is 
misleading  since  an  artificially  high  stiffness  was  used  at  this  location 
because  it  was  at  the  intersection  of  two  zones  of  the  smeared  pile  stiffnes¬ 
ses  . 

120.  When  shrinkage  was  added  to  the  analysis  with  ABAQUS,  the  loads 
were  redistributed.  There  was  a  s1 ight  increase  in  loads  under  the  chamber 
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area,  a  significant  increase  in  loads  under  the  wall,  and  a  decrease  in  loads 
at  the  outer  edge.  The  maximum  pile  load  occurred  under  the  wall  with  a  value 
of  360  kips/pile.  However,  the  previous  location  of  maximum  pile  load  still 
carried  a  load  of  350  kips/pile.  Therefore,  other  locations  can  be  considered 
critical . 

Cracking  potential 

121.  Comparisons  of  results  were  also  made  to  evaluate  the  potential 
for  cracking  of  concrete  in  the  L-17  analyses  reported  herein.  Modulus  of 
rupture  test  results  cited  in  paragraph  90  were  again  used  as  a  simple  cracking 
threshold  for  tensile  stress.  As  a  review,  the  modulus  rupture  values  used 
were  124  psi,  280  psi,  320  psi,  and  464  psi  at  1- ,  3-,  7-,  and  28-days  age, 
respectively.  The  following  conclusions  are  made  regarding  cracking  potential 
in  monolith  L-17. 

122.  First,  maximum  principal  stresses  computed  in  WES2DT  analyses 
which  included  gravity  and  thermal  loading  with  creep  reached  225  psi  in  the 
bottom  of  the  base  slab  (lift  1) .  These  stresses  developed  over  a  period  of 
more  than  21  days  as  a  reaction  to  thermal  expansion  of  lifts  2  through  5. 
Consequently,  the  peak  value  was  not  reach  until  the  concrete  in  lift  1  was 
more  than  21  days  old  and  had  attained  a  modulus  of  rupture  of  about  425  psi. 
Thus  the  maximum  tensile  stress  was  less  than  55  percent  of  tensile  stress 
capacity . 

123.  Secondly,  maximum  principal  tensile  stresses  computed  in  ABAQUS 
analyses  using  UMAT1,  which  employed  the  identical  modulus  versus  time  rela¬ 
tionship  as  WES2DT  analyses,  reached  300  psi  in  the  top  of  the  base  slab  by  5 
days  after  placement  of  lift  7  and  250  psi  in  several  locations  in  the  wall 
(Figure  48).  The  300  psi  level  occurred  in  concrete  that  was  at  least  15  days 


old  and  the  250  psi  values  developed  around  5  days  after  placement.  Even 
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of  rupture  attained  at  the  concrete  age  when  maximum  tensile  stresses  occurred 
124.  Finally,  maximum  principal  tensile  stresses  computed  in  ABAQUS 
analyses  using  UMAT2 ,  which  contained  an  excessively  high  modulus  function  for 
the  first  day  or  so  of  age,  reached  300  psi  in  the  bottom  of  the  base  slab  by 
21  days  after  placement  of  lift  1.  Stresses  were  only  about  100  psi  in  the 
same  location  and  time  when  creep  was  applied.  Stresses  in  the  wall  without 
creep  reached  around  200  psi  and  with  creep  were  about  100  psi  (Figures  50  and 
52).  Considering  that  the  computed  early-age  stresses  were  probably  high  due 
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Figure  42a.  Displaced  structure  from  gravity  turn-on  analysis 
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Figure  43b.  Displaced  structure  5  days  after  placement  of  lift  7,  gravity 
loading  only,  using  ABAQUS . 
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Figure  43c.  Displaced  structure  5  days  after  placement  of  lift  16,  gravity 
loading  only,  using  ABAQUS . 
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igure  UUa.  Maximum  principal  stress  contours  5  days  after  placement  o 
lift  4,  gravity  loading  only,  using  ABAQUS . 
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Figure  44b.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  7,  gravity  loading  only,  using  ABAQUS . 
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Figure  47b.  Displaced  structure  5  days  after  placement  of  lift  7,  gravity 
and  temperature  loading,  using  ABAQUS  (UMAT1) . 
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igure  48a.  Maximum  principal  stress  contours  5  days  after  placement  of 

lift  4,  gravity  and  temperature  loading,  using  ABAQUS  (UMAT1) 
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Figure  48b.  Maximum  principal  stress  contours  5  days  after  placement  of 

lift  7,  gravity  and  temperature  loading,  using  ABAQUS  (UMAT1) 
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Figure  49a.  Displaced  structure  5  days  after  placement  of  lift  4,  gravity 
and  temperature  loading,  using  ABAQUS  (UMAT2). 
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Figure  49b.  Displaced  structure  5  days  after  placement  of  lift  7,  gravity 
and  temperature  loading,  using  ABAQUS  (UMAT2). 
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Figure  50a.  Maximum  principal  stress  contours  5  days  after  placement  of 

lift  4,  gravity  and  temperature  loading,  using  ABAQUS  (UMAT2) 
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Figure  51a.  Displaced  structure  5  days  after  placement  of  lift  4,  gravity 
and  temperature  loading  with  creep,  using  ABAQUS  (UMAT2) . 
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Figure  51b.  Displaced  structure  5  days  after  placement  of  lift  7 
and  temperature  loading  with  creep,  using  AEAQUS  (I'M 
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Principal  stress  contours  of  structure  20  days  after  lift  13  is  placed,  gravity  and  temperature 
loading  including  creep,  using  WES2DT. 
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Figure  55a.  Displaced  structure  5  days  after  placement  of  lift  4,  gravity 
and  temperature  loading  with  creep  and  shrinkage,  using  ABAQUS  (UMAT2). 
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Figure  55c.  Displaced  structure  5  days  after  placement  of  lift  16,  gravity 
and  temperature  loading  with  creep  and  shrinkage,  using  ABAQUS  (UMAT2) . 
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Figure  56a.  Maximum  principal  stress  contours  5  days  after  placement  of 
lift  4,  gravity  and  temperature  loading  with  creep  and  shrinkage,  using 
ABAQUS  (UMAT2) . 
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Figure  56b.  Maximum  principal  stress  contours  5  davs  after  placement  of 
lift  7,  gravity  and  temperature  loading  with  creep  and  shrinkage,  using 
ABAQUS  (UMAT2). 


PART  XI:  CONCLUSIONS  AND  RECOMMENDATIONS 


125.  Throughout  the  course  of  this  study,  the  ABAQUS  computer  program 
seemed  to  adequately  capture  the  behavior  of  the  lock  monoliths  during  constru¬ 
ction.  Options  existed  in  the  program  to  adequately  model  the  phenomena  that 
occurs  during  the  incremental  construction  process.  While  at  the  present  time 
there  is  a  need  to  refine  the  aging  creep  model  used  during  this  study,  it  does 
appear  that  an  adequate  two-dimensional  incremental  construction  thermal  stress 
analysis  can  be  accomplished  using  this  computer  program. 

126  Accuracy  and  precision  of  results  in  finite  element  thermal  stress 
analysis  of  mass  concrete  structures  are  strongly  dependent  on  the  accuracy  and 
precision  of  material  properties  test  results  such  as  moduli,  specific  creep, 
shrinkage,  etc. ,  which  are  used  as  input  to  the  numerical  model  in  the  computer 
program  For  the  aging  creep  numerical  model  used  in  ABAQUS  in  this  study,  it 
is  c1  ar  that  the  early-age  E  versus  time  values  are  higher  than  those  values 
for  he  Locx  and  Dam  2 6  1  P  ■  concrete  being  modeled  This  led  to  calculation  of 
excessive  thermal  stresses  and  creep  relaxation  at  early  time  A  better  repre¬ 
sent.!'  ion  f  r  the  modulus  should  be  developed  and  verified  for  concrete  similar 
to  that  used  in  Lock  and  Dam  ?►>  P. 
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necessary  to  perform  an  incremental  construction  analysis,  preferably  with  an 
aging  creep  model,  in  some  instances,  if  practical. 

128.  Incremental  construction  models  better  simulate  actual  field  con¬ 
struction  conditions  than  simple  gravity  turn-on  analyses  in  producing  realis¬ 
tic  cons  true t ion  -  re lated  stresses  It  was  found  that  the  interaction  of  the 
manv  mechanisms  at  work  during  construction  can  produce  results  that  are  not 
obtainable  bv  simple  analysis  methods  More  research  is  required  in  this  area 
to  delineate  the  extent  to  which  incremental  construction  analyses  should  be 
required  in  structural  design.  It  has  become  evident  that  the  effects  of  the 
construction  process  should  be  considered  in  the  design  process,  because 
cracking  in  mass  concrete  structures  during  the  construction  phase  rr.av  alter 
design  considerations 

1 2  '<  As  pointed  out  in  paragraph  -*5,  the  costs  of  performing  ‘he  three 
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three-dimensional  analyses  of  incremental  construction  problems.  This 
formulation  was  developed  in  a  modern,  general  purpose  heat  transfer  and 
structural  analysis  code.  Every  effort  was  made  to  ensure  that  the  analysis 
concept  was  effective,  rational,  and  consistent.  With  this  requirement  in  mind 
it  soon  became  apparent  that  effective  modeling  of  key  parameters  affecting 
final  stresses  in  incrementally  constructed  mass  concrete  structures  were  pushed 
far  beyond  the  current  state  of  the  art.  These  key  parameters  include:  early 
age  (1  day,  2  day,  ...)  properties  such  as  shrinkage,  creep,  material  aging, 
and  cracking  strength.  Also,  the  concept  of  the  incremental  formulation  itself 
presents  some  unusual,  subtle  problems. 

132.  Currently,  research  programs  are  underway  to  better  define  and 
verify  these  parameters  through  carefully  planned  theoretical,  computational, 
and  experimental  research  programs.  The  intent  of  this  research  is  to  provide 
effective  and  efficient  design  and  analysis  guidance  to  the  field  engineers  in 
the  area  of  mass  concrete  construction.  The  approach  taken  here  is  a 
cooperative  research  effort  among  research  and  design  staffs  at  District  and 
Division  Offices,  major  universities,  private  engineering  firms,  and  Corp’s 


laboratories 
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APPENDIX  E:  ABAQUS  JOB  CONTROL  LANGUAGE  FILES 
USED  WITH  CYBERNET  SYSTEM 

FOR  TEMPERATURE  AND  STRESS  ANALYSES  OF  MONOLITH  L-17 


CYBERNET  JOB  CONTROL  LANGUAGE  FILE  USED  IN  ABAOUS  THERMAL  ANALYSES  (1,-17) 


/JOB 

ABAQ , T20000 , P2 . 

/USER 

/CHARGE 

PURGE, L7TT8/NA. 

PURGE, L7TPLOT/NA. 

PURGE , L7TREST/NA . 

DEFINE , TAPE12-L7TREST . 

DEFINE, TAPE8-L7TT8. 

GET , TAPE15-L7TN . 

GET , TAPE16-L7TE . 

BEGIN, , ABAQUS , I-L7TD,TEXT-N,APLOT-Y,USUB-Y, INSUB-LDFLUX . 
DEFINE , L7TPLOT . 

REWIND, NPFILEA. 

COPYBF , NPFILEA , L7TPLOT . 

EXIT. 

DEFINE , L7TPLOT . 

REWIND, NPFILEA. 

COPYBF , NPFILEA , L7TPLOT . 

/EOR 

/EOF 


L7TN-node  file 
L7TE-element  file 
L7TD-conunand  file 
LDFLUX-heat  subr . 


CYBERNET  JOB  CONTROL  LANGUAGE  FILE  USED  IN  ABAOUS  STRESS  ANALYSES  (L-17) 


/JOB 

ABAQ , T30000 , P2 . 

/USER 

/CHARGE 

ATTACH , TAPE17-L7TT8 . 

PURGE, L7SRES/NA. 

PURG  E , L7  S  PLOT/NA . 

DEFINE , TAPE12-L7SRES . 

GET , TAPE15-L7SN . 

GET , TAPE16-L7SE . 

BEGIN, .ABAQUS , I-L7SD , TEXT-N , APLOT-Y , USUB-Y , INSUB-L7SUBS . 
DEFINE, L7SPLOT. 

REWIND, NPFILEA. 

COPYBF, NPFILEA, L7SPLOT. 

REWIND, TAPE35. 

COPY.TAPE35 . 

REPLACE , TAPE35-L7ST35 . 

EXIT. 

DEFINE. L7S PLOT. 

REWIND, NPFILEA. 

COPYBF , NPFILEA , L7SPLOT . 

REWIND, TAPE3 5. 

COPY , TAPE35 . 


L7SN-node  file 
L7SE-element  file 
L7SD-command  file 
L7SUBS-modulus  subr. 


E2 


APPENDIX  F:  ABAQUS  HEAT  GENERATION  SUBROUTINE  DFLUX, 
FILE  "LDFLUX"  USED  IN  2-D,  L-17  ANALYSIS 


ooooooooo 


SUBROUTINE  DFLUX( FLUX , TEMP , KSTEP , KINC , TIME , NOEL , NPT , COORDS , 

&  JLTYP) 

DIMENSION  COORDS (3) ,Q(23) ,T(23) 

COMMON  /ELDEF/  STIME(316) 

DATA  ENTIME/960 . 1/ 

DATA  NQ/23/ 

DATA  T/  6.,  12.,  18.,  24.,  30.,  36.,  42.,  48., 

&  60.,  72.,  84.,  96.,  120.,  168.,  240.,  288., 

&  360.,  408.,  480.,  528.,  600.,  720.,  960./ 

DATA  Q/0. 017385,  0.015158,  0.011682,  0.010309,  0.00800, 

&  0.C06862,  0.005795,  0.004910,  0.003263,  0.002897, 

&  0.002211,  0.002059,  0.001357,  0.000980,  0.000628, 

&  0.000507,  0.000368,  0.000305,  0.000241,  0.000214, 

&  0.000201,  0.000175,  0.000118/ 

DATA  STIME/ 

&  28*0 . 0 , 14*0 . 0 , 28*120 . 0 , 28*240 . 0 ,28*360 . 0 , 22*480 . 0 , 16*600 . 0 , 

&  16*720.0,16*840.0,16*960.0,16*1080.0,16*1200.0,16*1320.0, 

&  16*1440.0,16*1560.0,16*1680.0,8*1800.0/ 

C  ********************************************************************* 
C  VERSION  OF  USER  SUBROUTINE  "DFLUX"  USED  FOR  MONOLITH  L-17,  2-D  MODEL 
C  ********************************************************************* 
C  VARIABLE  DEFINITIONS¬ 
ENTIME  -  END  OF  RELATIVE  HEAT  GENERATION  TIME  +  SMALL  TOLERANCE  (HR) 
NQ  -  NO.  OF  HEAT  GENERATION  RATE  POINTS 
T  -  RELATIVE  HEAT  GENERATION  TIME  POINTS 

Q  -  HEAT  GENERATION  POINT 

C  STIME  -  VECTOR  CONTAINING  PLACEMENT  TIM  .  FOR  EACH  ELEMENT 
C 

C  FLUX  -  HEAT  GENERATION  RATE  RETURNED  TO  PROGRAM 
C 

TREL  -  TIME  -  STIME(NOEL) 

IF (  TREL. GT. 0.0. AND. TREL. LT.ENTIME  )  GO  TO  10 
FLUX  -  0.0 
RETURN 
C 

10  CONTINUE 
FLUX  -  0.0 
DO  20  I-l.NQ 
J  -  I 

I F (  TREL.LE.T(I)  )  GO  TO  30 
20  CONTINUE 
C 

WRITE(6 ,35)  KSTEP, KINC, TIME, NOEL 

35  FORMAT(/ , "  WARNING  -  PASSED  THROUGH  DFLUX  WITHOUT  ASSIGNING", 

&  /,"  FLUX.  STEP  -",15,"  INC  -",15, 

&  /,"  TIME  -" , F12 . 2 , "  ELEMENT  -",15) 

RETURN 

30  FLUX  -  Q(J) 

RETURN 

END 


APPENDIX  G:  ABAQUS  USER  SUBROUTINE  MPC, 
2-D  MULTIPLE  POINT  CONSTRAINT 
VERSION  USED  WITH  MONOLITH  L-17 


o  n 


nv 


SUBROUTINE  MPC (UE , A , JDOF , N , JTYPE , X , U , NMPCE) 


DIMENSION  A(N) ,JDOF(N) ,X(6,N) ,U(6,N) 

COMMON/COUNT/KINC , MINC , KITER , MITER , FATIME , ATIME , DATIME , 

1  CTIME , DCTIME , DTIME , DDTIME . HTIME , DHTIME , DDTPRE , DATPRE , HTIM1 , 

2  DHTIM2 , EXFAC , KSTEP , KCUTS , MCUTS , NUMBER , LSHAF1 , LCUTBK , DTNEWS , 

3  KITGEN , KMDINC , TTIME , DTMIN , DTMAX , MITEIG , MITXXX , STIME , DSTIME 
4 , TPREV , TNEW , TOLD , TEND 


C 

C 

C 

C 

C 

C 

C 

C 


********************************************************** 

SUBROUTINE  MPC  USED  WITH  L-17  2-D  ANALYSIS.  MPC  IS  MERGED  WITH 
SUBROUTINE  UMAT1  OR  UMAT2  IN  FILES  "L7SUBS"  AND  "L7SUBN,  RESPECTIVELY 
********************************************************** 


FIX  NODES  ACROSS  TOP  OF  VOID  DURING  PLACEMENT  OF  LIFT  8 
STEPS  29,30,31,32,33,  AND  34 


IF(KSTEP.EQ. 29. AND. JTYPE. EQ. 101)  GO  TO  10 
IF(KSTEP . EQ . 30 .AND .JTYPE . EQ . 101)  GO  TO  10 
IF(KSTEP.EQ. 31. AND. JTYPE. EQ. 101)  GO  TO  10 
IF(KSTEP.EQ. 32. AND. JTYPE. EQ. 101)  GO  TO  10 
IF(KSTEP.EQ. 33. AND. JTYPE. EQ. 101)  GO  TO  10 
IF(KSTEP . EQ . 34 .AND .JTYPE . EQ . 101)  GO  TO  10 


10 


NMPCE  -  C 

RETURN 

CONTINUE 


B1 

B2 

B3 


X(l,3) 
X(1 , 3) 
X(l,l) 


X<1, 2) 
X(l,l) 
X(1 , 2) 


UE  -  (B2*U(2 , 2)  -  B3*U(2,3))/B1 


A(l) 

A(2) 

A(3) 


-  B1 

-  -B2 

-  -B3 
JDOF(l)  -  2 
JDOF ( 2 )  -  2 
JDOF(3)  -  2 


RETURN 

END 


ABAQUS  USER  SUBROUTINE  MPC,  2-D  MULTIPLE  POINT  CONSTRAINT  USED  WITH  L-17 
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APPENDIX  H:  MATERIAL  USER  SUBROUTINE  UMAT, 
VERSION  "UMAT1" ,  2-D  MODULUS  ROUTINE  WITHOUT  CREEP 
USED  WITH  MONOLITH  L-17 


APPENDIX  I:  ABAQUS  USER  SUBROUTINE  UMAT, 
VERSION  "UMAT2" ,  AGING  CREEP  MODEL  WITH  SHRINKAGE 
USED  WITH  MONOLITH  L-17 
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NPLAST-  ARUS.2)  AA-  TH/57.29578 

NPL1-  NPLAST / 1 00  C«  COS(AA) 

NPL2-  (NPLAST-NPL1*100)/I0  S-  SIN(AA) 

NPLJ-  NPLAST -NPL1 • 100- NPL2*1Q  10  CONTINUE 


SOIOtN)-  SIG1  tPHl  l'DEPS  XK-XKlDELTM'O  25 

50  CONTINUE  ‘0  CONTINUE 

i- 1  1(1)-  2  0*SG(6tN)!l  0  XKN-  1  .  0/XK 

SI J<  2 1 »  -  SC  (  6 1 N  )  /  }  0  CALL  COE E I XKN , T , AF  ACT , TEMFAC) 

SIJ(5)«  SC(fcfN)IJ  0  CON-CHEEP / AFACT 


